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Theorems concerning the properties of Kennard packets have been extended over all 
the solutions of the wave equation of a harmonic oscillator. Generalizing these theorems 
to the quantum theory of fields, non-spreading wave packets with clasiscal properties 
have heen obtained. 


Introduction 


The outward spreading of wave packets associated with a free particle, makes its 
realistic interpretation impossible. However, it was known for a long time past that 
in the case of a harmonic oscillator a non-spreading wave-packet may be obtained 
(Kennard 1927); the factor ensuring, so to say, the stability of the packet is the harmonic 
potential. The „Kennard packet“ which is obtained by introducing the time variable 
into the general Gauss function presents very interesting features. Its centre vibrates 
in conformity with the laws of classical mechanics, and its mean quadratic deviation 
in position (measuring at the same time the width of the packet) is a periodic function 
of time. The oscillator is, however, only a very special problem in the framework of 
non-relativistic quantum mechanics; it may also seem that the properties of the Kennard 
packet are intimately connected with the *Gaussian^ form of the initial packet; it is 
not clear why just such a form of the initial packet is distinguished in theory from 
all the others. These considerations do not allow us to connect with the Kennard 
packets any deep seated physical meaning. ; 

On the other hand, however, it is known that such fundamental physical systems 
as quantized free fields (in the relativistic theory) may be considered as consisting of 
harmonic oscillators in the sense of non-relativistic quantum mechanics. It may be 
expected, therefore, that in the case of these fields non-spreading wave packets with 
classical properties might be found. The costruction of such packets is one of the 
chief objects of the present paper. 


* A short account of фе main formulae of this paper was published іп the Bull. Acad. Polon. 
Sci. (1954). 
(24155 
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Before entering upon the subject of the “field“ packets, we shall give a prelim- 
inary discussion of the mathematical mechanism assuring the non-spreading of 
the Kennard packets. In the first part of our work we shall consider some properties 
of general time-dependent solution of the harmonic oscillator wave equation, and give 
some theorems generalizing in a certain sense the properties of Kennard packets over 
all the solutions of the wave equation. By specializing these theorems, we shall get an 
infinite family of packets (connected with the stationary state functions) with 
properties quite analogous to “Gaussian Kennard packets“. Afterwards the theorems 
applying to one harmonic oscillator will be extended over quantized fields. 

We shall make use of some properties of unitary operators given in the paper 


by Infeld and Plebaüski (1955): “On a Certain Class of Unitary Operators“?. 


§ 1. Notations and auxiliary formulae 


The Dirac symbolism (Dirac 1947) will be adhered to throughout the paper. 
Let x, p, w be the position, momentum, and circular frequency of a harmonic oscillator 
of mass m. We shall employ the dimensionless quantities ; 


E, = (m w|h)* x, & = (moh)-^ р. (1.1) 
df df 
In quantum theory &,, £, are thus dimensionless real operators, which, due to 
[p, x] = — hi, satisfy the commutation rule 
[$i 62] = $1 6. رغ غ س‎ = (1.2) 
df 


The Hamiltonian of the oscillator H = 1/2 [p?/m + mo? x?] expressed in terms 
d 
of Ё is Н = oh Ху, where X,— 1/2 (22, + é). Thus, if instead of the time 2 we make 
df 


use of the dimensionless quantity t = wt, then, in virtue of 
df 
ee. 

H|¥>=ih>|¥> (1.3a) 
the vector describing the state of the oscillator will be given in terms of an arbitrary 
initial state | p> as | 

Roto ri^ | m | (1.3b) 

We shall denote by {п >}, n = 0, 1, 2,... the orthonormal system of eigenkets 
of the operator X,: 


Xo In > = (n + 12) п> (14a) <п'|п" > = 0, (1.45) 
и | Лаа . (1.4) 
п=0 


We may choose the phases factors of the kets [n so that their representatives obtain- 
ed by means of the eigenbras of the operators Ё. (EJE, > = EX lé» = 


| 1 This paper will be hereafter referred to аз IP. 
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= 6 (Fı — £j) ) аге given by the symmetric formulae 


«|n = U, (&) e" (1.5) 

<& | п> =U, (&) e "4 (1.5Ь) 

U, (2) = n-^ (27n1)-* Н (2) e (1.5с) 
df 


In these formulae H, is the n-th Hermite polynomial and the consistency of 
(1.5a) with (1.5b) follows from the well-known formula 


+00 
1 > 
= ~iba = [п 
әл С Г» (b) = GA (a) 


From the properties of the generating function of the Hermite polynomials it is 
easily seen that one can also obtain the functions U, (z) defined in (1.5c) by means 
of the (symmetric) generating function: 


exp [— 1/2 22 + V2 zu — 1/2 w?| = x^ y и" U, (z)] Үл!. (1.6) 
п-0 


It follows from (1.4с) and (1.4a) that the operator 675% appearing іп (1.3b) can 
be expressed as 
со 
e iX. У) |n > е "+ < n | VIS 414) 
п-0 
and hence, in accordance with (1.5a—b), the representatives < & | e 4244635 > 
of this operator are 


e Le KIr &. А (1.82) 
< & [e| & > = K[r — zl, 8, 6] e"*, (1.8b) 
< & e^ | ê > = Кт + ae, £j, ёе, (1.8c) 

ЕЕ, Fal, (1.84) 
where 
со 4 à 
Kir, гш] У, (г) U, (ш) et a (1.93) 
df n=0 É 
This last function may also be written as? 
К [mz, z, w] = і" 6 [z — (—1)" w), еі (1.95) 
{же 


(sina) ^ expi/2 [(z? + w?) ctga — 2(—1)" zw cosa] 


Үл (1.9с) 


K[mz + а, 2, w] tm 


where 
Cie fas I т = 0, +1, + 2,... 


y 


2 The way of deriving (1.9b) and (1.9c) from the definition (1.92) is given in IP. 
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82. The periodic time-dependence of the probability amplitudes 


According to (1.3b), the probability amplitudes in the position and momentum 


spaces: 0 (т, £) = <&,|W>, Р (t, £) = <, |У are expressible in terms of 
af df 


their initial values О, (£j) = < £& | o >, Po (£j) = < &, | Po > as 
df df 


199 + со ; ^ p. га 
О (т, 6) = [< Ele" 2 > 42, Q (6) = |< E E > 48, Р/6) (2718) 


+оо +оо | 1 қ қ 
P(x, 6) = | « &le Lo а РАВ = f < Ве Е > а, (Е) (225) 
From these formulae it follows easily, thanks to (1.8), that 
OE = Р( 22, & ее, (2.2а) 
PT EN OE eee jest (2.2Ь) 


These formulae show that the probability amplitudes in the position and momen- 
tum spaces stand in an analogous relation with each other as the position and the 
momentum of a classical oscillator. In classical theory, using the dimensionless ,,posi- 
tion“ &, and „momentum“ £, defined in (1.1), we get 2, from £, by bringing it a quarter 
of a period back. 

The same analogy applies to the mean values. We shall denote the expectation 


value of any operator 4 by А, = <Y | A | ¥ >. Let (2) be an arbitrary function of z. 
df 
In virtue of (2.2a), we have 


FE) = f 4810 (r, &) PFE) = f 4&1Р(т—л{2, E f (ê) = Ев (23) 


In particular: 
5, = dee 52, Ж Bart 
and hence an analogous relation between the mean quadratic deviations? follows 
4 5, =A BUS 


Apart from this quarter-of-a-period difference in phase between the quanti- 
ties in the position and momentum spaces, it is evident from (1,93), (1.8), and (2.1) 
that the functions Q(t, £j), P(t, £j) are periodic in time. 

Thus with an arbitrary initial condition | Y, >, a uniform spreading of the initial 
packet is out of question. 


a uu AS 


! Ab, = [< E SLE P14, 


ДЕ = [<¥ | (& — V | E | WS) РІЗ. 
* 


MEL. Ма 
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The mechanism of the periodicity of the functions Q(z, Е) and P(t, &,) makes 
it possible to further elucidate the nature of the relation (1.9b); making use of this 
relation in (2.1), we find according to (1.8) that 


О (kr, £) = i^ Q [(— 1)* £j] (24a) 
Q[(k+ У] л, &) = itt Р, [(— 1)* Е] (2.4b) 
P (kr, £) = i-* P, [(—.1)* &] (2.4c) 
P ]) — 35) л, &] = i^^ Q, [(— 1(* E (2.44) 


Thus we see that probability amplitudes oscillate, so to say, between their initial 
values (related by the Fourier transformation). This is clearly shown in the follo- 
wing table. 


О (т, $1) = i% Py (ë1) 1 
P (t, ё) = 20 (22) г? Po ($2) 
83. Some general properties of the oscillator wave-packets ' 
Let us now consider the following unitary operator 
ghi ای کے‎ ой е- کے‎ АЙ. о. ЫЙА, (3.1) 


where é, are arbitrary numbers‘. From (3.1) it is easily seen that the representatives 
of the state | Wy > = её — Êê) | V, > are given by 
df 


Og: (E) estt pau Oy (EE) exp i [88h E2 
df | 
D (E) = = s EE, = Py (& — &) exp i[ — 6, +6, 6/2, (3-Б) 


where О, E = < &4 >, Р(&) = = xh | Wo > are the representatives of | $ >. 

Let us yii. assume ұла we know OTE) Gele IS, 
Р (т, &,) = < £, |e *° | Wo >. Then the probability amplitudes corresponding to 
the initial state | V, > may easily be found in consequence of the Heisenberg equa- 
tions for an oscillator. We have 


e E, eX — cos т E, + sin v £y ебх, e "№ = — sin т ё + cos тё. (3.38) 


4 Тһе proof of the correctness of the identity signs, as well as of formulae (3.2), may be found 
in ІР, 89, point e. 
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It follows that 


e i ей ZI ei X — gii) (3.3b) 
where the numbers n are given by 
& = Е cos т + 2 sin т, p = -- i sin T + F cos T, (3.3c) 
so that they are solutions of the classical equations of motion of an oscillator 4. & = 
E 55% Zi = & = — 8) with arbitra initial conditions. 


б [са with (3.3b), we have therefore 


| >= ети | ЧЛ > = e i die 5) | % = 
df 4 


gi 756) e i" | Po > = ghi 58) | v — (3.3d) 


Thus finally, in virtue of (3.3d) and (3.2), we find that the representatives of the 
state | W > associated with the initial state | Wọ > are 


о.) = <ê 
P(E HIDE =Plr, & – Е,) ехрі[– 65,4 5, 6,12], (3.45) 
df : j 


where Ез E are the classical quantities from (3.3c) depending on т. 

The relations (3.4a—b) generalize in a certain sense a property of the Kennard 
packets — the motion of its center according to the laws of classical mechanics — 
over all the solutions of the equations of motion of an oscillator. On co uparing (3.2) 
with (3.4), we see that the introduction of the time-dependence produces a »revival« 


in the classical sense of the parameters £j, &, which displace the initial packet with 
respect to the centre of attraction of the oscillator. Thus the classical quantities 


& (т) characterize in a certain sense the position of the centre of the packet in 
phase space. The quantity ЖР & Ë appearing in (3.4) may be also interpreted physically 


by saying that w = 3F (5,2, - zi E has the meaning of the classical action connected 
d 


with the oscillator’. 


5 Let H (pj, 4; t) be a second-degree homogeneous function of the variables р; and q;. Because of 


te ta 
the dle idewtity FP Ji 3 do ни = Se Í | purto Жаргы dr 
df A i 27% 2 Эр; 2 94; 


: ӘН ӘН 
On using the Hamilton equations: dp; = — ehe dt, dq; — PUN dt, we get W = 5 үн | P;dq 
қ і і 


; 1 t 
+ 1/2 q;p; — 1/2р; in| ut S Pidi | ? This auxiliary theorem is being used in the text. Evi- 
- 1 

к> 1 
. dently w is dimensionless. 


ELSE exp HUE S BLEU р Дд: 


Шы с а. 
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A posteriori the relation (3.4) may be derived in a simpler manner. Assuming 
that О (т, £j) satisfies the Schrödinger equation, one can confirm by substituting 
О (т, &) into this equation from (3.4a) that Q is also its solution; it may be also di- 
rectly seen that Q verifies the initial condition (3.2a). 

| One could also obtain P (г, £) from О (е, 6) by means of a Fourier transforma- 
tion. : 


Let us now consider the unitary operator 


Е, (a) = e , X, — 1 (5 & + & Ê) (3.59) 
S df ағ 2 
This operator acts іп a quite definite manner on the state vectors?; namely, if the 
functions 0, (£i) = < £y >, Py (&) = < Е, | V, > are given, then after putting 
| Vo > = Ra (а) | Yo >, we get 


Q(&) = < $1 | Po > =a" Qo (£a?) 
P(E) = < &lyo = a* Py (& a). 

Let us now assume that we know the probability amplitudes О (т, £j) = 
< & AUR id Wy >, P(t, é) = < Ё, | e "X | цуу > connected with the initial state 
EZ >; than the probability amplitudes corresponding to the initial state 
| Yo > = R, (а) | V, > follow immediately. Indeed, the operator identities 


е "К. (а) = eih e log b (т) As [b (т)] е їагсіва "tg TX, (3.6a) 


е-і R, (a) = eth log C (т) В, [с (т)] e- arts ate tX, (3.6b) 
where | 
b (т) = (a? cos? т + a^? sin? т), c(t) = (a^? cos? т + a? sin? т)% (3.6c) 
df gf 4 
then hold. They сап be verified as follows: differentiate the left hand sides and the 
right-hand sides with respect to тапа transpose afterwards the operators to the right 
taking account of their way of acting on the £a's (formulae (3.3a) and (3.5b)); thus 
the equality of the derivatives of both sides of (3.6a —b) may be proved; the validity 
of the identities in question for т = 0 is self-evident. 
In accordance with (3.6a—b), we find on using (3.5c— d) 
Eod 
ift, 1980 (т) 
(3.7a) 
и, E log C (® 


(3.7b) 


О (t, £) = < & | e ** |B > = b7 (v) Q [arctg a~? tg t, £i b? (z)] е 


P(t, &) = < &e I | Po > =с (т) P [arctg a? tg T, бс 1(т)] е 


where b (т) and c (т) are defined by (6.3c). 


6 See ІР, $9, point 6. 
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A posteriori (3.7) could be proved by verifying that Q (т, &;) is a solution of 
Schródinger's equation when О (т, £j) is its solution and P (т, 2) from (3.7b) is the 
Fourier transform of О (т, £j) from (3.7a). 

The relations (3.7a —b) show that a »change of scale« in the initial condition (the 
parameter а in (3.5c— d)) leads concerning the time-dependent wave function to 
a somewhat different time-dependence, and to time-dependence of the change-of- 
scale parameters. 

Bringing together the theorems expressed by (3.4) and (3.7), we come to the 
following result: Let Q (т, &) = < &|e 7 | Po >, P (т, £) = < Eafe fep 
be known probability amplitudes in the position and momentum spaces corresponding 
to the initial state | yy >. Then to the initial state 


е 343 о Е 1 
| Wo > = eie bo етйова > (В+) | и, > (3.8) 


there correspond the probability amplitudes 


Qs E) = « & Lec |, > = 57A (9) О [aretg ават, ш aul T 


exp i E & (т) — w(t) + 1, (61 — & (9)! = log b (т) + & ы» | (3.9a) 


LC аа у СЕ Б atig, $ 4 х 


exp i | - 25 Е, (т) + w (T) + Yo (2 — È, 0) < - log ере à ба | (3.9b) 


where 


& = & "cos т + 5, sin т, & = — É віпт + E сов т (3.9c) 
df df 


are solutions of the classical equations of motion of the oscillator with the initial con- 
ditions &,, &,, the quantity 

w (1) = EN (т) Ea (T) — È Èa) (3.9d) 
being the classical action for d motion and b (т), c (т) being given by 


b Ня ка cos? т + a^? sin? т)”, c (т) = (a^? cos? т + a?sin-? т)? (3.9e) 
df 


Equations S may be directly verified in a similar manner as it was done for (3.7). 
The states | Y > = е7 | V, > whose representatives are given by (3.93— b) 
are obtained from | Y^, > by successive application of unitary operators; therefore, 
if < Yo | > = 1 (which we assume from now on), < ¥ | ¥ > is also equal to unity. 


It might be added that on putting 4, = = А, А, == SALA urs. 


we get 


“DS 
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on =b (т) а агсіва” ?tgc A à (т), Se C (т) Eo arctga*tgz -} ^ (т), (3.10а) 
and that the analogous mean square deviations in | Y >> are connected with the devia- 
tions in | V > by the relations 


ДЕ (v) = b (т) А Ё, (arctg a^? іст) 
A & (т) = c (т) Л é, (arctg a? іст) (3.10b) 
Thus we see that our general wave-packet (3.9а—Ъ) has in a certain sense the 
properties of a Kennard packet. The centre of the packet moves according to clas- 
sical laws. The change-of-scale parameters (in comparison to the initial packet corre- 


sponding to | Wo >, b (т) and с (7)) are periodic functions of time and are simply 
connected with the quadratic deviations. 


§ 4. Special wave-packets 


Particularly interesting is the application of the general formulae (3.9) to the case 
of a state | W > (associated with О (т,&1) and P (т, &,)) depending on the time т only 
through a phase factor. This property is obviously shared only by the stationary 
states | V, > = e X | п > = e *@+%) | n >. From (3.9) it is easily seen thanks 


df тати 
to (1.4) and (1.5) that if the numbers A £,,, A 6, are interconnected by the relation 
AE, A E, En M, NOD 
the probability amplitudes in the position and momentum spaces of the states 
| > = gr y ысы) [л B, (n + Ye) |а o (4.1b) 


df 
are respectively aom lE S>, р, = EE huis: 

Aeg (О ap iss e] Ез = 
nV» = d^ E 2 == ИЖС 
qalt, &1) я = E E,}/ VE exp i | &#& —и 

— (n + ¥f,).arct A £y t ul (гек. ty log ДД (т) — & Ê la + n ml 
А Ла 5 1 ә 1 I n Qu в 4 
` (4.22) 


A A [45 0 ]- ue. Аё (1) er ЧЕ CMS 
Pn (т, £3) EM al £3 % Шш A p 2 


P 


- (n+ мав щт + Û (E Б о CO) +È J2 — na | 
к (4.2Ь) 


The classical quantities Е, A are here defined by (3.9c); in virtue of (3.10a), 
we find that in the case under consideration these quantities represent mean values 


of the position and momentum in the considered states; Е апа Б have thus 
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the meaning of initial values of these average quantities, w (т) defined by (3.9) is the 
classical action. Because of (3.10b), the quantities 


А &, Kies AG £i? cos? т + A 52 sin? т)”, ЛЕ, (t) = (A 12 віп? t + A £j cos? т)“ 

(4.2c) 
may be identified with the mean square deviations in the considered states. Thus 
A Evy Ace have the meaning of initial values of these deviations; because of (4.1a) 
only one of these quantities may be chosen arbitrary. Thus, the physical significance 
of all the quantities appearing in (4.2a—b) have been explained. 

Thus we arrived at a family of wave- -packets characterized by a discrete para- 
meter п and continuous parameters ct Ей 4. £2; these packets generalize in a natural 
way the »Gaussian« Kennard packet; (it is obvious that A £i, can be arbitrary given 
for every n). 

It may be easily shown that our packets exhaust all the possible solutions of the 
oscillator wave-equations (in the position and momentum spaces) whose time-depend- 
ence lies explicitly (i. e., neither by means of the mean quadratic deviation nor the 
expectation values of position and momentum) in the phase factor alone. 

It may be also worth mentioning that after fixing the parameters EAR. а = E. 
(п +1/)~% for a whole family of our wave-packets, the functions (q, (т, &1}}, 
(p, (т, &,)} form obviously for every moment v complete orthonormal sets in the 
position and momentum spaces. 

To facilitate the discussion of the properties of the wave packets (4.2a —b), 
we introduce the function 


US EY zx efe- Шет, 71 (4.32) 


We assert now that for every n 
lim 6, (x — x, Ax) = à (x — x), (4.3b) 
4x20 
‘so that the parameter Ax (denoting the mean square deviation in x when д, is con- 
sidered as the probability density in x) measures, so to say, the width of the 6, func- 
tion. 


_ Indeed, because of (1.5c), lim ô, (x — x, Ах) =0 for x zx. We have also 
Ax +0 


+00 +00 


ilm f die òn (x — % 29/04 = lim J di! U(x!) f [=+ 5 А = f) 
and (4.3b) follows. 
We shall denote the probability densities in the position and momentum spaces 


connected with | V, > by —, (&, т), 5, (E; т). They are given respectively by 
Ta @ьт) = ê, [fı — &, Ain (J, с, (En т) = 9, [Ea — & A &, (0] (4.4) 
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These relations show that the centres of our packets move classically (ғ, = 4, (z)) 
and that the parameter which must vanish in order to make the function д, coincide 
with the Dirac function varies periodically with time — thus an indefinite spreading of 
our packets is impossible; in the course of time they undergo only periodical 
change in shape. Because of the Heisenberg relations (cf. (4.1а)) it is obviously impos- 
sible so to choose the parameters as to make either 7, 0:7, coalesce with the Dirac 
function. For example, if we take a very small initial А s, then LES will be very 
large, so that the initially very sharp packet 4, will have, according to (4.2c), a width 
pulsating in time between a small AEs, and a large /4&,. 


For the state | W, > the product of 9 mean square deviations in positions and 
momentum is 


A Ge (т) 4 Ды. (т) = [(n + 3/9? + (42 — A £22 sin? т cos? т)|% (4.5а) 


and hence it depends оп time if only A &, +A Рб 
The smallest deviation appears for т = Ёл |g, and the largest one for т = 


= (k + 15) л}: 
A Fig A Fog | ening = + Ya A bin A Fon | ve (ықы, Ma (4 En 4) (455) 


It is also easily seen that the expectation value of the energy in a state | W, > is 
< Xo > = "h (Ef + ED (LES + AES) 
It is composed of a classical part and a »fluctuating« part. At fixed 8, апа z the 
„fluctuating“ member will have, in virture of (4.la), a minimum for En = № +: +4, 
EL ч e 
In this case: (1) The deviations in position and momentum for|V,- do not 


depend on time; (2) the shape of the packets for the probability densities is also time- 
independent, only the centers of the packets oscillate in a classical fashion; (3) < x, > 
= n + 1, + 1» (& + 2), hence the expectation value of the energy is the sum of 
the energy of the n-th stationary state and the classical energy. In this case the proba- 
bility amplitudes are simply. 


an (n £) = 0, (& — ор [E & — w — (n 4-39 v — & 2 + пл] (4.78) 
p. (ту &) = U, (& — &) expi[— Ею — (n + Ya) 7 + 5,12 — пл (АЛЬ) 


One gets the smallest value of the energy, or the < Ху > associated with such 
a state when 6. Е approach the limit zero; 4, and p, from (4.7a—b) change then into 
probability aniplitudes related to the n-th stationary state, as expected. 

Let us now return to the general case. Hitherto we have been working with dimen- 
sionless quantities. However, if in (4.4) we restore the initial x, p and t, we obtain ex- 
pressions for the probability densities which, due to (4.3b), possess the following 
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properties: 
lim g, (x , t) = ê [x — xy(t)], lim „(руд = ô [p' - рыд)! (4.8) 
A0 =o 
where 
хы (t) = xg cos wt + po[mo * sin wt, ры (t) = — max, біп wt + cos wt 
df 


with (mo/h)” E = xy and (moh)^ s. = ру. The numbers 4 éin and A £, are 
arbitrary but fixed, and subject only to the limitation (4.1a). 

It is to be emphasized that the limiting process (4.8) is correct independently of 
the value of A £j, / A E» It is also interesting to note that the right-hand side ot 
(4.8), i. e. the classical probability densities, do not depend on n. Therefore it would 
be wrong to suppose that on the quantum level only one wave packet corresponds to 
each given classical motion ,,(¢), р„(@) of the oscillator. In the limiting case h — 0 all our 
wave-packets (4.2a —b) turn into a classical motion, which is not the case, however, 
for quite arbitrary wave-packets О (т, £j), P (v, £j). The property lim o (х, 0) = 

%->0 


= Û (x — x (t)) holds only for our wave-packets associated with stationary states. 
To summarize we can say that our wave-packets (4.2a — b) generalize the Kennard 
wave-packet in a similar sense as the 6, (x — x, Ax) functions generalize the Gaussian 
distribution. Their properties are very similar to those of a classical material point in 
a harmonic field of forces. They are marked out from all other solutions of the Schró- 
dinger equation by the exceptional admissibility of the limiting process lim o (x , 275% = 
ћ > 0 
= 0 (x — х, (t)). Also the fact is important that for all n-s at fixed values of the 
parameters т, &, &, а = A &, (п + 1/2)” * our packets form complete orthonormal 
sets of functions in the position and momentum spaces’, 


85. The wave-packets in the case of field theory 


The harmonic oscillator and the related wave-packets are of utmost impor- 
tance not so much in the case of an oscillator considered from the point of view of 
classical mechanics (a rather academic example) as in the application to relativistic 
quantum field theory. | 

It is well known that if a general field y4(x, t) satisfying a linear relativistic equa- 
tion (the superscript A reforring to the mode of transformation of the field) is expanded 


? These facts suggest the idea of looking for an extension and improvement of the usual statistical 
interpretation, so as to elucidate on the one hand relations between our wave-packets and classical mecha- 
nics, as well as the relations between the д, (х — x, A <) functions and the calculus of probability, and on 
the other hand the minimum properties of the stationary states. 

We shall not dwell upon these questions here. The aim of this paper is only the construction of 
generalized Kennard packets in accordance with the usual interpretation and above all the discussion 
of the significance of wave packets in the framework of the field theory, 
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within a periodicity box into a Fourier series, i. e. into plane waves, then the Hamilto- 
nian of the field reads 


ШЕ МО EC COO X HA (5.1) 


The quantities P, and Q, stand in linear relations with the Fourier coefficients of 
$^(x,t) and their canonical conjugate momenta z^(x,t). The index n describes 
here the degrees of freedom (we include in it not all the different degrees of freedom 
connected with the A’s, but only the effective degrees; e.g., in the case of the elec- 
tromagnetic field P, and Q, are connected only with the transverse vibrations of 
the field). The classical equations of motion can be then put in Hamiltonian form, 
as follows 

ӘН ӘН 
20; N^ Era 


In quantum theory P, and Q, have to be replaced by operators satisfying the 


n 


(РА == (5.2) 


usual commutation rules and the equation for the vector representing а state can 
be written in the form 


2 
ЗЕТ (5.3) 


The ket referring to а stationary state of the field тау be represented in terms of the 
eigenkets of the operators H,, which we shall denote by | N,> (Н, | N, > = (М, + 
ERN as: — j 


Шу = МЫ, exp [eios (МЫ (5.4) 

The kets | №, > have representatives км іп the Q and P spaces as 
«Q, IN, = агу, Uy, [04а] i5; a, = ODÊ (5.52) 
«P, | N, > = fz” Una (РУВД і В, = (ho) (5.5) 


From the above relations and the completeness of the set of kets |N; ,... N,...— = 
А, | N,> it follows that the representatives of the ket | "n ЖО O 


can. E obtained from its initial value V, (01..., Q,...) as: 


W (Qi... no fag. fd. - K [t5Q y+ 0,23 1. Onee] 001, Qn) 


where d 
ORTA roe ЫЕ О. > өнүн ОН с ай) г з= (8:06) 
df мМ 
Paying due regard to the ME of X and П її is easily seen that in (5.6b) 
one can rearrange these орсо, as follows 


Él Qu QU] =D < Q, | N, > em NO, (5.60) 
Ў df n М, 


-0 
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By (5.5a) and (1.9a) we see that K may be expressed by means of the previously 
investigated K function as 


K [t5 Q,...5 0*..] = Па K[o, t Ola 0/0, (5.6d) 


It follows immediately from the above formula that all the previously discussed 
properties of the wave functions of a linear oscillator can be taken over into the field 
theory. | 

Among others, the following generalization of formulae (3.4a-b) are true®: If 
V (Qj... Q,...; t) satisfies the initial condition ¥ | = Yo (Q1... Q;....), then to the 


г=0 
initial condition 


, , о , о 1 . ست‎ м 2; 
Py (Qi... О» ...) = Po 101 — 0ь-->0.- Om Texp m X (Q, P, CAPIO 
4 ; (5-72) 
where Р, 0, are arbitrary numbers, there corresponds a state with the representative 


0002, ..., Quy 12) = р101 — Qu ы On — Qu ird өр У (03Р, — 15,0) 

" (5.7b) 

Here 0,0), Р (1) are solutions of the classical equations (5.2) with the initial conditions 
On Р п 

Obviously an analogous formula exists іп the P, space. Тһе matter is a little differ- 

ent with the generalization of formulae (3.7). This generalization succeeds only in 


the physically rather uninteresting case when all o, are for all п equal to one and the 
same constant. | 


On using (5.7a-b) in the case of a |W> corresponding to the state 
[N,, ... М„...; t >, we get the wave packet 


V ess (5 Оль -.. Qu) = Haz ^ Us, (0,- O,())/a,] 


x exp y У (0,2,0 — FPO 0,0 — һи, (s. 4 i t+ Na лу) (58) 


and hence it is clear that for the expectation values in the state Bom Eros 


we always have, irrespectively of the values of N,,... Np ... : 


<0,>--0,(0 =н Рр) (5.9) 


Since the field operators z^(x), ф“(х) depend linearly on P, and Q, it follows 


8 For simplicity we omit the operational formulation of the theorem. 
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from the above that the expectation values of the field operators for АМ г> 
are identical with the classical fields (to which in the language of Fourier transforms 
correspond the initial conditions P,, Q,). This is true for arbitrary combinations of 
М, ... N,.... Concerning the energy, the case is a little different. The expectation value 
of H in the considered state is 


<H> = > һ о (М, + 15) + 1, > (Р? + oig) (5.10) 


Thus the energy consists of two parts, the first one representing the energy of 
the stationary state described by |N,... N,..., t >and the second one being in an evi- 
dent manner identical with the classical field energy. 

In the above considerations the solutions have not been restricted to particular 
fields; they are valid therefore for all fields whose energy can be put in the form (5.1), 
i.e., in all normal cases. Of course, our work was of quite a formal nature, as we did 
not pay attention to any convergence questions, and thus from a mathematical point 
of view many objections could be raised (concerning, e.g., the convergence of the 
infinite product (5.64)). 

Formulae (5.7a-b) make sense when the field is contained in a periodicity box 
and we have to do consequently with an enumerable number of degrees of freedom. 
It may be presumed, however, that the contents of these formulae remain valid quite 
generally, independently of the enclosure of the field in a periodicity box. 

It is well known that one can circumvent the difficulties inherent in the use of 
the periodicity box by working in a representation in which the field functions y“(z) 
are c-numbers and the state vector is a functional of these functions. We shall now 
show the validity of the content of (5.7a-b) in the above representation, considering 
for concreteness the case of a scalar meson field. 

Тһе Schródinger equation for the vector representing a state of the scalar field reads: 


| 2 | 
[esto + (rot te 95 [р> =й 5, > (5.11) 
where the operators 2(x), p(x) satisfy the commutation rule 
Іл(%), p= —ihàs G^ — x" 


Let ||ІФІ> be the common eigenket of the operators ф(х) (for all x), i.e., 
9 (х) |[p']>= 9 G) | [p] 2 ; the eigenvalues ф' (x) of the operator ф (x) are c-number 
functions of 7. Their ensemble characterizes the eigenket of the operator g (x), 
depending „functionally“ on the function g' (x). Consequently in the ket symbol 
ф'(х) has been enclosed in square brackets. 
We define now the functional of state as 
V [p] @) $ «[v] | v7 (5.13) 


reserving once for all square brackets for denoting functional dependence. 
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The commutation rule (5.12) may be satisfied by identifying л(х) with the oper- 
ator of the functional derivative:? 


л (x) = —thod/6d o (x) 


In the new representation equation (5.11) becomes 
д? , , , . д и 
Ul ЕЕ — hrc? $o Gy + (o)? + и? ф | Pip (0] = ih = V [p] (t). (5.14) 


We shall now prove the following theorem: Let У [y] (t) be an arbitrary fixed solution 
of equation (5.14); then 


7 [o"] (t) F V [p' — Ф (t)] (0) ° exp ИЕ p a(t) — "o (5.15a) 


where о (x, t) л (x,t) are arbitrary solutions of the classical equation of motion 


дф д 
= т. AE pus 
еа , (4 — и?) p = 2: ^ (5.15b) 


and // (t) is the corresponding classical action of the field (is also a solution of equation 
(5.14)). 

The easiest way of carrying through the proof is by substituting (5.15a) into 
(5.14). We shall go into some details to illustrate the methods of computation by means 
of functional derivatives. We have: 


6% | 
Meri ames 


| 20197 — 900] 99 RD, ‚ 99 


‘co 


aW ; E 5.16 
= 24 ' exp a [...] V [p — v] (9. MEUM 


In the first member-on the right ~~ 20 : ) denotes differentiation of W with respect 


to the time whereby only the explicit (not by means of g(t)) dependence on time is taken 


‚ 9 Ф1 ® Оф’ + en'l] — 2 [g'] 


ОШ) (x) d E ТЕЙ (x) зх 
Re(x) 
together with its derivatives outside the region К, (2) and on its surface. a is an inner point of the 
region R,(z) and the distance of any two points of this region is always smaller than =. We define the linear 
operator “615 iz) Ы) the relation 


‚ Where 7)(х) is an arbitrary function vanishing 


е es [ord ааа <> 
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into consideration. We get the second term on the right by applying 2 Q [o (4)] = 


p(x, 4) д 
dx e ni DP The symbol бо) (%) denotes that the functional derivative 


um to be к with respect to ф' and then, after differentiation, g' has to be replaced 


by o' — ^. The origin of the third term is obvious. 
poe it is easily seen that 
ô? i 62 T ES: 
ӛр” (3)? ; V [v'] ] (0 = vg y [бсо] dy? (2) 2 
2@{8Ч]р—Ф](@) — i д jl LO ir ob pt i 
z [e] dax 9' t - "-о|()---- -- 
а aou) Ft ETO LM) ы exe C] 
i ó? Te S 
= exp > eas or" Y [p' — p] (t) (5.16b) 
2i д 
+ 5 л 55 У | t А EL V [p' — 9] (t 
"rp 3 [y — A 


Thus when Фр" (4) is substituted into (5.14), the left-hand side of that equation 
can be written, thanks to (5.16b), as follows 


oe 52% vo 
= са даа. * оқа 272 
HW = exp >= ^ aT JE L| drl- 2 209 (ауу 2hic? x 3 (с (3) + сл 


+ (rey + и? ез) TP [p — v] (t) (5.16с) 
This expression may be further transformed into: 
r^ А Pe | , E 1 FR 
HY = exp +. 4 e Бәз? pr tne - gts - $] (t) 


52 бек TRIN E ف‎ 2992802 - 2772 
көрігі [4=| 2 ihc “гү? o o (9) +202 фф 0% 


52% x У [4 — 91(0 
According to the assumption us the first set of terms on the right-hand side 
of (5.164) may be replaced by ih —— = 0 V [p' — g] (t) exp x [...] Taking this fact 
into account, we see by comparing (5.16a) with (5.16d) that in order that 
ih 2 y [91 (¢) = H ¥ [$] (t) be valid independently of the shape of the functional 
W [p'] (t) the following relations must hold 
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Ф - oe ôF [р — $] (t 
— ih E ——ihchm (өзгені of Eme (5.16e) 
dn Wl s Lm 
- fap gy + 0 ے‎ DH yo + u?* Ф 9) 
+ = fa x G л? — (үф)? — и? Р (coefficients of ¥ [p — $] (t) (5.16f) 


These relations may be also written in the form 


S 
+ - ШЕ А гаа 70 | 
0 oo 


This relation must be identically true for any function being an arbitrary argument 
of the functional. Hence, it will hold only in the case of equations (5.15b) being satis- 


nrw ^ s п 
fied and if W (t) = z IE fos (c л — (Yq)? u2 о?) +e const, аз well aS 


о (dg/dn) | г. = 0. Thus  (t) represents in fact the classical action, and in order 
for equations (5.14) to be satisfied by V [g’] (t) the classical equations (5.15b) must 
be verified by лапа о (relation (5.16b) is identical with the first of equations (5.15b)). 
The condition рф (dp/dr) |; may be satisfied in two ways: either the class of the 
functions Ф’, on which our functionals may depend, is a class of functions vanishing 
at infinity — then ф would be a quite arbitrary solution of (Г] — и?) ф = 0 — ог, 
when nothing can be said of the functions ’, we must postulate (dg/dn) |. оо = 0. In 
such a manner, however, plane waves would become inadmissible. 

The above-given calculations prove thus the theorem contained in (5.15). It 


seems rather obvious that if the state ¥ (¢) > to which V lv] (t) corresponds is sta- 


tionary, the expectation value of Ф (x) and л (x) for the state Y (t) > will be identical 
with y (x, t) and л (z, t) irrespective of the stationary state in question. (This intuitive 
mode of reasoning is justified in a certain sense by the relations (5.9), as this property 
ought to be preserved in the limiting case of an indefinitely large periodicity box). 

The generalization of the theorem contained in (5.15) over the cases of other 
fields does not present any serious difficulties. In general, one may expect (because 
of the splitting of any quantized free field into harmonic oscillators) that, if the field © 
p(t) (where A denotes tensor or spinor indices) complies with linear relativistic 
equations, the following theorem will be true: 


On Certain Wave-Packets 293 


If the state functional У [ ф’4] (t) satisfies the wave equation of the given field, 
the same equation will hold for the functional 


V (410) = Y [p^ — 9^ (010) exp т IY | d ха HA (1) — 70| (5.17) 
А о 


where p (x,t), л“(х,) are arbitrary solutions of the classical field equations, and W (2) 
is the classical action connected with them. 

The author wishes to express his thanks to Professor L. Infeld and Professor 
ХУ. Rubinowicz for their continuous interest and encouragement. 
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SELF-DEPOLARIZATION AND DECAY 
OF PHOTOLUMINESCENCE OF SOLUTIONS 


By A. JABLONSKI 
Physics Department, Nicholas Copernicus University, Torun 
(Received August 4, 1954) 


Self-depolarization of photoluminescence of solid solutions {depolarization by migra- 
tion of the excitation energy from one luminescent molecule to other luminescent molecules 
of the same kind) is treated in a similar way as the quenching of photoluminescence of 
solutions was treated in a previous paper. The problem is, however, even more compli- 
cated than that of quenching, and thus some additional simplifying assumptions have to 
be made. The model of a luminescent centre adopted here is the simplest one. The excited 
centre is assumed to consist of an excited luminescent molecule surrounded by an “active 
sphere“, which may contain one or more unexcited luminescent molecules of the same 
kind. If there are unexcited molecules within the active sphere, the excitation energy can 
be transferred to them and migrate several times from one molecule to another. The 
probability of an energy transfer is assumed to be proportional to time and be the same 
for transfer to any of the unexcited molecules belonging to the same centre. It is assumed 
to be nil for transfer to unexcited molecules outside the centre. 

АП the centres of a system (i. е. in our case, of a solid solution) are divided in groups 
according to the number of unexcited (depolarizing) molecules belonging to a centre — 
a centre is ascribed to the group A if it consists of one excited and &—1 unexcited mole- 
cules. Expressions are derived for the probability as a function of time of finding the 
excitation energy with a molecule primarily excited by absorption of light. These expres- 
sions are different for different groups of centres. Hence, the process of depolarization 
runs differently in different groups. There is e. g. no depolarization at all in the group 
k= l; 

Expressions are obtained: (i) for the decay of the particular components of photo- 
luminescence, and (ii) for the degree of depolarization (or that of polarization) as a func- 
tion of the concentration of luminescent molecules. The latter expression is compared 
with an experimental curve published by Chauchois. 


81. Introduction 


Some processes produced in molecules by interaction with tlieir neighbours, 
as quenching or depelarization of photoluminescence, occur in the condensed state 
in a way differing from that in which similar phenomena occur in gases. In the case 


(295) 
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of gases the probability of quenching, or that of depolarization, is simply proportional 
to the number of effective collisions per second and is the same for all excited mole- 
cules of the same kind in the same macroscopic conditions. The situation is different 
in the case of the condensed state, say, in the case of liquid or solid luminescent solu- 
tions, in which the probability of the above processes is in general not the same for 
all luminescent molecules but depends on the configuration of molecules surrounding 
a particular excited molecule. In order to describe adequately the phenomena occuring 
in this last case, one has to divide all the excited centres! of a solution into groups 
according to the configuration of interacting molecules (quenchers or depolarizing 
molecules) and to ascribe to centres of each group a particular value of the probability 
of a process being caused by interaction. The behaviour of each group of centres 
has to be described by a separate equation. Such was the way in which the quenching 
phenomena were treated in a previous paper (Jabloüski 1954). In the present paper 
the self-depolarization of photoluminescence and its influence on the decay of partic- 
ular components of photoluminescence have to be treated in a similar way. Only the 
case of solid solutions will be considered below. Even this case, although much simpler 
than that of liquid solutions, is sufficiently involved as to require some simplifying 
assumptions to be made. Some results of the paper on the theory of quenching (Jablon- 
ski 1954) have to be used here. 


82. Migration of excitation energy between molecules of the same kind 


The migration of excitation energy in an isotropic luminescent solution from 
a luminescent molecule primarily excited by absorption of light to other molecules 
of the same kind is considered as the cause of the so called “self-depolarization“ 
of photoluminescence (cf. Jabłoński 1936; Vavilov 1942, 1943; Pringsheim 1949, 
page 378). The molecules excited by energy transfer have practically an isotropic 
distribution of orientations, and hence the photoluminescence emitted by them is 
practically completely depolarized. 

The probability of an energy transfer depends of course on the configuration 
of the unexcited luminescent molecules surrounding the excited one, and thus is 
` different for different centres. 

Let us first consider the following problem: Two identical molecules 4 and B in 
a fixed distance from each other are given. Let one of them, say molecule A, be excited 
at the moment ¢ = 0. What is the probability of finding the excitation energy at 
molecule 4 at t > 0, if the probability of the energy transfer from one molecule to 
the other during dt is assumed to be udt, u being constant for a fixed distance between 


——— —— HOUR 


1 ,Photoluminescent centre“ means here a luminescent molecule together with neighbouring 
molecules surrounding it. 
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the interacting molecules 4 and В? The excitation energy can wander from A to B 
and back from B to A several times during t. If the number of transfers is nil or even, 
the excitation energy will be found at molecule A, if this number is odd the molecule 
B will be found excited. It can be shown that the probability W (Г) that i transfers 
take place during £ is given by? 


Р()- 82 eit, (1) 
Hence, the probability W, that at the moment t the energy may be found with 


the molecule A (i.e. the probability that the number of transfers which occured during 
t is even ог nil) is 


z 2m —2 


The probability of finding the excitation energy with B is 


ZR peti inh pees pm (3) 
= -- o-Ht = د‎ 
А РЕ. Я (Оле)! =e “sin h ut 2 . 


For 2 > co, W, and W, > 1/,, as should be expected’. 

Eqs (2) and (3) cannot be applied directly to real solutions, since there are in such 
solutions centres with very different configurations of the unexcited molecules, to 
which the excitation energy can be transferred, whereas the above equations apply 


only to pairs of molecules with a quite definite mutual distance. 


2 This expression is given here without proof, since, as my colleagues mathematicians informed me, 
such a proof can be found in many treatises on the theory of probability. I am sorry to be unable to 
quote any of these treatises. 


3 It is perhaps worth noting that the mean values W, and JV, of W, and М, which would be оБвету- 
ed if the molecule A would be excited several times, or if we had to do with a system consisting of 
several identical pairs of molecules 4 and B, are given by 


1 —2ut 1 1 
= elem D dcl mus 
2 2(1-- ux) 


and 


со 
ШЫ ondes e ыы чап 1 
тет a уа д” 
0 


where т denotes the mean life time of an isolated excited molecule. 
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The situation in real solutions is very complicated indeed. Thus it seems to be 
advisable to introduce into the theory some simplifying assumptions. Accordingly, 
the model of a luminescent centre adopted here is the simplest possible. The excited 
centre is assumed to consist of an excited luminescent molecule surrounded by an 
„active sphere“; the excitation energy can migrate only between those luminescent 
molecules which are within the active sphere; the probability of transfer of energy during 
dt from one molecule to another one is assumed to be the same for all pairs which can 
be formed from luminescent molecules belonging to the centre and equal to udt; 
the probability of transfer of energy from a molecule inside the sphere to a molecule 
outside it is assumed to be nil. A more refined model, the ,,shell model* introduced 
in previous paper (Jabłoński 1954), should correspond better to real conditions, 
it would lead, however, to expressions too complicated to be of any practical use. 

Let us consider a centre with k luminescent molecules, from which one, say 
molecule A, be excited at t = 0. We look'for the probability W_,(t) of finding the 
excitation energy with A at t > 0. W(t) must satisfy the following conditions: 
Я (t) = 1 for t = 0, and W ,(t) ^l/k for t > oo (the probability in question must 
become the same for all & molecules belonging to the centre, and thus must become 
equal to l/k for t > co): 

We now assume 


dW 1 
Ka а (ma - 1]. (4) 


where a is a constant to be specified later. The solution satisfying the conditions just 
mentioned is 


1 1 
У. 1-1-2) EUR (5) 
Since W4 = 1 for t = 0, Eq. (4) can be written as 
k—1 
аўд == - Т а. (6) 


On the other hand, the probability that any transfer of excitation energy occurs 
in a centre containing one excited and Ё —1 unexcited molecules is (k — 1) times 
that of a transfer to a single unexcited molecule. Thus for t = 0 (W = 1) 


аў == — (k — 1) ий. (7) 
From Eqs. (6) and (7) there results 
а = ku (8) 
Putting this value into Eq. (5), we obtain 


1F (Fa) er 


ro = (9) 


lied. 
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For k = 2 Eq. (9) becomes identical with Eq. (2), the last equation thus being 
a particular case of Eq. (9). Since Eq. (2) seems to be well founded, the identity of 
both equations for Ё = 2 may be considered as justification of the assumptions from 
which Eq. (9) is derived. 

The probability of finding the excitation energy with one of the molecules B 
(i.e. with one of the primarily unexcited molecules) is 


(10) 


We shall make use of Eq. (9) in the next sections. 


83. Some general remarks on polarization of photoluminescence 


Let an isotropic photoluminescent substance (body) be excited by plane-polarized 
light vibrating parallel to the z — axis of an rectangular co-ordinate system x, y, z, 
and let L,, L, and L, be the relative intensities of the components of the emitted pho- 
toluminescence along these axes. The symmetry of the problem requires that L, = L}. 
We put L, = L, and Г, = Г, = L, , and normalize the intensities so that 


ГЕР, = 1. (11) 


If the emitted photoluminescence is polarized, then L,,= L}. 
The degree (rate) of depolarization of photoluminescence is defined as 


Ly 
==, (12) 
Ly 
and that of polarization P 
Ны ніл de 02s (122) 


Lye by 1+ 0 
Denoting 


Li -L =r T (13) 


and taking Eq. (11) into account, we may write 


VEM (14) 


and 


L = (15) 


From (12a), (14) and (15) we have 


(16) 
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and 
2I 
= 17 
eae (17) 
From (12), (14) and (15) we obtain 
l—r 
XE КЕ 18 
аз; d 
and 
1-0 
== 19 
А 20-1” 09) 


Let a luminescent solution be suddenly excited at і = 0. If there are depolarizing 
factors in the solution, the momentary values of |r|, | 1 — | and | P | will decrease 
with time. In usual experimental conditions when a luminescent substance is illu- 
minated with the exciting light continuously, or when the duration of illumination 
is long enough, mean values of the above quantities are observed (averaged over 
the time £ extending from 0 to 2o). These mean values can be calculated for particular 
cases if the mean life-time of the excited state of the molecules in question and the 
variation of r (or o, or P) with time are known (cf. also Jablonski 1935 and 1936). 


84. Decay of polarized photoluminescence of solid solutions 


As already mentioned, the problem of self-depolarization is even more involved 
than that of quenching. An additional complication arises e.g. from the fact that the 
migration of the excitation energy from one molecule to another leads in general to 
a change of configuration of the centre, and thus to its transition to another group 
(one of the exceptions being the case of centres constituted of two luminescent molecules 
only; the migration of the excitation energy leaves then unchanged the only parameter 
characterizing the configuration, i.e. the mutual distance of the molecules). If the 
transitions from one group of centres to another would be taken into account, the 
behaviour of the system had to be described by a system of interdependent differential 
equations, in which case the calculations would become very tedious. 

As already mentioned (cf. § 2), we adopt here the simplest model of a centre, 
that of an active sphere surrounding the molecule in question. In order to describe 
the behaviour of a luminescent solution, i.e. of a system of luminescent centres, we 
have to divide all the centres of this system into groups (cf. 8 1). We divide them as 
follows: group of centres with one luminescent molecule (the excited one) only (group 
1), group of centres with pairs of molecules, one of them excited (group 2), and so on; 
generally, group Ё is assumed to consist of k — 1 unexcited and one excited molecule 
(k = 1, 2, ... oo). The orientation of the excited molecules is not taken as a criterion 
of belonging of a centre to a certain group. Thus, the migration of the excitation energy, 
although causing depolarization of photoluminescence, is not assumed to produce 
transitions of centres from one group to another. 
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The behaviour of each group of centres has to be described by two equations: 
an equation describing the variation with time of the population of the group, and 
another describing the process of depolarization. The equations of the first kind were 
given and discussed in the paper on quenching (Jabłoński 1954). For the case of solid 
solutions excited suddenly at the moment t = 0 they read: 


d 
Te mtn Wy ps. (1 1,22. ос) (20) 


with the solutions 


Py = py, exp L.-(y + Wy + Wy, (Е = 1, 2,... >) (21) 


Here у is the probability of a transition leading to the emission of photoluminescence, 
W,— the probability of inner quenching plus that of quenching by foreign mole- 
cules (which is assumed for simplicity to be the same for centres of all the groups), 
7, — the quenching probability by unexcited luminescent molecules‘ i.e. by the same 
molecules which cause the self-depolarization, and p, denotes the population of the 
group Ё normalized in such a way as to fulfil the condition 


со 
У рь г. 1, 
Е-і 


where p,, are the values of p, at t = 0; p,, is thus the probability of finding a centre 
with (2-1) unexcited molecules at ¢ = 0. Let the volume of each centre be v cm? 
and the concentration of the unexcited molecules (which is almost always practically 
the same as the total concentration of all luminescent molecules) be n cm 3. Putting 
vn = v, we have (cf. Smoluchowski 1907) 


ны у : dp 22 
Pko (k 1) ! е ( ) 
From (21) and (22) we obtain (cf. Jablonski 1954) 


yk-1 


Р гулбу SEP [—(y + W; + Wa t], (fos со) (23) 


Eq. (23) describes the decay of population of each group of centres in the case of 
sudden (very short) excitations at 2 = 0. 
We have now to consider the equations describing the process of depolarization 
` of photoluminescence by migration. | 
Let us divide each group of centres into two sub-groups: sub-group А of centres 
in which the excitation energy is with the molecule originally excited by absorption 
of the primary (polarized) light, and the sub-group B of centres in which the excitation 


ғ; 


4 This kind of quenching is often called ,self-que nching“. 
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energy is with one of the primarily unexcited molecules. The photoluminescence 
emitted by sub-group 4 is in general polarized. Its degree of depolarization, the so 
called „fundamental depolarization^^ may be characterized by the corresponding 
value ry of г (cf. Eqs. (13), (18), and (19)). 

If the probability of excitation energy transfer to an unexcited molecule does 
not depend upon its orientation (as it is assumed here), the distribution of orientations 
of the excited molecules belonging to the centres of sub-group B must be uniform, 
and thus the photoluminescence emitted by it must be completely depolarized‘. 
Thus, only centres of sub-group А contribute to the value of r in each group. According 
to Eq. (9) the probability that a centre of group Ё belongs to the sub-group 4 is 


1 + (k — ење 
а. 


The variation of г with time for the group Ё is given by 


1+ (k — l)e-*^ 
rh = ro Ға) = А (24) 
If Eq. (19) is taken into account, (24) can be written as 
1—0 1 + (k— 1) e 
Wer cu OE crt ages (25) 


— 1+ 20, k 


Eqs. (24) апа (25) show that, for finite k, г, does not vanish when t - оо, i.e. the 
depolarization never becomes complete. This property distinguishes the depolari- 
zation by migration from that by other mechanisms, such as Brownian rotation of 
the luminescent molecules in liquid solutions or depolarizing collisions in gases. 

The relative intensities of the components L,, and L,, of photoluminescence 


emitted at various times t by the centres of group k and normalized so as to fulfil the 
condition LZ, + 2 L,, = 1 аге (cf. Eqs (11), (14), (15), and (24)) 


1 + 2r 1 1 р — Nemh 
Гы mre iren est k ps | (14а) 
апа 
КЕ 1 1--(Е-І1)е-ен 
я (1-6 ет. (15а) 


5 In most papers оп the polarization of photoluminescence the notion of „degree (rate) of polari- 
zation" is rather used than that of ,degree (rate) of depolarization". The last notión seems, however, 
to be preferable in some cases and was used in several papers (Jablonski 1935, 1935a, 1936) 

$ Vavilov (1952) assuming a specific mechanism of energy transfer arrives at the result that the 


migration probability depends on the mutual orientation of the molecules. Presumably this effect, if 
. . LP : 
taken into account, would not change significantly the results of this paper. 
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If Len and Г, are multiplied by the relative population p, of group k, there result 
the relative intensities of the photoluminescence components ры = p, E, and p,, 
= p, Lp, with respect to the intensity of the photoluminescence emitted by the whole 
system: 


Ры = Pk Ly (26) 
SO 200. oon l — @ 1 + (k—1) exp (— kut) 
a ea sah кы 
and 
Pki = Pr Lra (27) 
Eg ii 1 — oo 1+ (#—1) exp (— ky t) 


The notations used in (26) and (27) are those previously used in (21), (23), (24) 
апа (25). It is well known that the „effective sphere“ for depolarization is much larger 
than that for quenching. It is to be noted, however that the notion of „active sphere“ 
used in this paper differs from the familiar notion of the „effective sphere“ in so far 
as the probability of a process with an interacting molecule within this sphere is here 
assumed to be finite, whereas such a probability should be infinite for the 
familiar ,,effective sphere“. We are thus allowed to admit the same volume of the 
active sphere both for depolarization and for quenching, provided we adequately 
assume different probabilities for these processes within the sphere. 

Тһе decay of the components of photoluminescence emitted by particular groups 
of centres proceeds according to (26) and (27) if the photoluminescence is suddenly 
excited at t — 0. 

The decay of the components of photoluminescence emitted by the whole system 
(by all groups) results from summing up the corresponding components emitted by 
all the particular groups: 


РИ РЬ | 

T (28) 
LE TU al {а Шафер be) —(у + Wy + Wa) 
i Е" 2 чаа ^^ 
апа 


puo дұры (29) 


k=1 


SENS cae: ‚Каре. -F W+ Wi) d. 
2 (mE k pigra A 
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From Eqs. (28) and (29) one obtains easily an expression for the decay of any other 
components: | 

Ра == Pj COS а + p, sin? а, (30) 
where a is the angle between the direction of the component in question and that 
of the electric vector of the primary light. 

As Eqs. (28), (29) and (30) show, various components decay in a different manner, 
which thus differs from that in which the total photoluminescence emitted by the 
whole system in all directions decays. Only the component making an angle а == 
54°44’ with that of the electric vector of the primary light decays according to the 
same law as the total photoluminescence: 

The laws of decay were already considered in two papers (Jabloüski 1935 and 
1936) for cases in which the probability of the process of depolarization is the same 
for all centres of the system (depolarization by Brownian rotation in liquids and depo- 
larization by collisions in gases). The expressions for the decay of p, and p, obtained 
there are simpler than those given by (28) and (29), р and р, being respectively simply 
equal to a sum and a difference of two exponential terms. It is worth noting that, 
whereas p [pj — 1 for t + oo in the cases considered in former papers, this is not 
so in the case of depolarization by migration considered here. 

Let us examine this question somewhat closer. The momentary degree of depo- 
larization, p(t) = p | [pj, is equal to the fundamental degree of depolarization o, at t = 
= 0, i.e. at the moment of the sudden excitation of the system. With increasing t 
о (t) gradually approaches its limiting value o (20), which, as already mentioned, 
equals 1 for the cases considered in the previous papers, but in general о (>) = 1 
in the present case. This limiting value can be easily calculated provided W,~ 0 
(no self-quenching). From eqs. (28) and (29) we obtain in this case 


nS UE ( 3: — 
E MUERE TA vacui GET TM 


Prete 31 

ЭО) a EE EE ee 

ee то Чоу»; IS a e PN 
АЛЕТ rcp 


which is equal to 00 for v = 0, and equal to 1 for v = co. 


85. Degree of polarization 


The degree of depolarization of the photoluminescence of an isotropic solution 
excited by continuous illumination is | 


[s dt 

а em (32) 
Гр 
0 
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For pj and p, as given by (28) and (29) this becomes 
oo 
те EE ER (1—00) (5-1) | 


£k! УИ, W, y+ 0 + Wet Еш 
м) 
Е sed 2 (1—00) (6-1) 
e k! ui. yt Wet ZEN 
If we put t = 1(у + Й,), this can be written as 
urea | 3ш ито [k(1 + 2ш) — (1—о)] | 
£2 (Е—1) 111+ Wet ku) vg (1 + Wa to) [1 + (Wr + ku) ту] 
Panter . (33а) 
А | 3 ите [k(1 + 200) + 2(1—0%)] | 
= (k—1) ! [1 + (Wr + ku) to (1+ Wa to) [1 + (We + ku) ту] 
The degree of polarization is given by 
(Ри = ра) dt 4 
120022 ХЕ CORR НИНИ а х, 
Гоп +p.) d 


0 
By means of the equalities 1 — oy = 2P)/ (1 + Po) and 1 + 200 = (3 — Р,) / (1 + P) 
one obtains from (33) and (34) 


1 6—1 
Ye 1) nl comm Pastore 


pecgpuu (85) 
yk-i Bis SPIER, P, (k—1) | 
= (kD! [y+ Fr + WF, yc + Wat іш 


Let us remark that the volume of the active sphere which has to be ascribed to given 
luminescent molecules is to a certain extent arbitrary. It depends, ceteris paribus, 
on the adopted value of the probability of migration и: the larger this value, the 
smaller the active volume (v = v/n, where п cm is the concentration of the lumines- 
cent molecules) must be assumed in order to fit the theoretical expressions to the 
experimental results. The above arbitrariness is a characteristic feature of the simpli- 
fied model of a centre used in the present theory. We may derive profit from this 


d 
arbitrariness and put u =y + Я, = — . We thus obtain instead of (35) 


со e [| 1 m k—1 | 
(5-11 Ё+1+ ьт 
peen sae ope ТЕ KR (35a) 
ә-і [3k — (k— 1) B, =I P, | 
м И И 
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If, besides, W,~ 0, the last equation becomes 


Ao pki Е-і1 
ene drum 1 


Р--3Р,-- c (36) 
У (E— на аса МЕ 


Remembering that у = vn, where v cm? is the volume of the centre, and n cm3 
the concentration of the (unexcited) luminescent molecules, we see that Eqs. (33), 
(33a), (35), and (36) represent the dependence of the degree of depolarization (or of 
that of polarization) on the concentration of luminescent molecules. 
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Degree of polarization of fluorescence of fluorescein sodium in glycerol versus concentration. O experi- 
mental points of Chauchois, A points calculated by means of eq. (36). 


The expression (36) in which the self-quenching is neglected (W, = 0) may be 
compared with the experimental results of Chauchois. Fig. 1 shows the curve published 
by Chauchois (degree of polarization of fluorescein sodium in glycerol versus concen- 


tration) together with points calculated by means of (36). Two of the calculated points 
had to be fitted to the experimental results. 


The agreement is quite satisfactory. 

In the general case the self-quenching cannot be neglected, and thus Eq. (33) 
or (35) has to be applied instead of (36). There are e.g. cases in which the degree 
of polarization increases with concentration. Such a peculiar behaviour is most probably 
due to the influence of self-quenching (Feofilov and Sveshnikov, 19407). 


It may be noted that if in (36) all terms but the first two are neglected, this equation 
can be written approximately as 


7 СЕ. also Levshin (1951). p. 274. Possibly the ,,self-quenching“ by polymerization plays an impor- 
tant réle in such cases. 
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1 = (A JR k 3m 2 | 
> | y 


The form of (37) is similar to the familiar form given by Perrin (1929) to his 
expression describing depolarization of fluorescence caused by Brownian rotations of 
luminescent molecules. (37) can be valid, if at all, at best for у <1. 


КРАТКОЕ СОДЕРЖАНИЕ 


A. Яблонский, Концентрационная деполяризация u затухание фотолюми- 
несценции твердых растворов. 

К теории концентрационной деполяризации фотолюминесценции твердых 
раотворов был применен тот же метод, который уже раньше был применен авто- 
ром к теории тушения фотолюминесценции. Получены выражения для затухания 
отдельных компонент фотолюминесценции и для зависимости степени поляри 
зации фотолюминесценции от концентрации раствора. 
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The problem of the Hamiltonian formalism in the case of first order equations of 
motion iseinvestigated on a dynamical model with a finite number of degrees of freedom. 
The possibility is shown of introducing canonically conjugate variables to the generalized 


coordinates g; by the usual definition pj = E, in spite of their dependence on 4» 
: q; 

provided that a suitable symmetrization is Е The relations between the 475 and 
p's are interpreted as a canonical transformation independent of time. The canonical 
quantization with the aid of the correspondence principle (exploiting the properties of 
the Poisson brackets) is proved to be possible in this case in spite of the mutual depen- 
dences of 4 and p. The correct canonical commutation relations with the factor 1); are 
obtained. The admissibility of treating the relations between the p's and q's as a canonical 
transformation is restricted to only one of theinfinitely many Lagrangians giving the same 
equations of motion, i. e. differing by terms of the form of a time derivative, namely to the 
one symmetrical with respect to integration by parts. This form of the Lagrange function 
(connected with the canonical commutation relations possessing a factor 1/;) is thus distin- 
guished. Moreover, two other special forms allowing an immediate elimination of one half 
of the variables q are distinguished too; then the remaining canonical variables, already 
independent, fulfil the traditional commutation relations without the factor 1/;. These 
considerations are generalized for the known cases of fields obeying first order field 
equations (Dirac and Kemmer fields). 


I. Introduction 


The problem of a dynamical system, with a finite or infinite number of degrees of 
freedom, described by first order Lagrange equations offers specific difficulties while 
performing a transition to the canonical formulation. As was shown by Burton and 
Touschek (1953) in two special cases (constructed within the framework of quantum 
mechanics of material points), the consistency of the Lagrange equations with the 
operator equation Q — ih [Q, H] is possible only when the MED MM 
between the generalized coordinates 4, and momenta р, defined as Tee 


(309) 
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are assumed in the form ([q,, р] = x ды, differing by the factor 1/ from the traditional 


form of the canonical commutation relations (established by the correspondence method 
of Dirac, or resulting, in the general case, from the variational principle of Schwinger 
(1951)). In order to overcome this difficulty and to deduce the canonical commutation 
relations Burton and Touschek used the variational method of Peierls (1952). However 
this metod is not indispensable. It was shown by Schwinger (1953a) on a general model 
(with a finite number of degrees of freedom) that this problem can also be solved 
by a suitable application of his variational principle; it must be taken into account 
that in the case of first order equations the p,'s are linear combinations of the 9,5 
and not independent, canonically conjugate variables. The proper, independent 
pairs of canonical variables are obtained by Schwinger after some transformations 
of the time dependent terms in the Lagrange function, giving the complete decompo- 
sition of the variables into two kinematically independent sets ,,of the first“ and ,,of 
the second kind“ (equivalent respectively to the division into the cases of Bose and 
Fermi statistics). The number of variables of each kind is shown to be even, and they 
can be grouped in pairs (possessing a character of canonically conjugate variables) 
for which the canonical commutation relations result from the variational principle. 
These considerations applying to systems with a finite number of degrees of freedom 
allow an immediate field-theoretical generalization where an analogous difficulty 
exists for the case of first order field equations. The general theory of these fields, 
based оп Schwinger’s variational principle is given in his ,,Theory of Quantized 
Fields II“ (1953b). 

It seems interesting, however, to investigate also the problem of systems obeying 
first order equations of motion from another point of view, based on Dirac’s correspond- 
ence method of canonical quantization, using the analogy between the ,,classical“ 
and the „quantum Poisson brackets“. Using, for simplicity, а model with a finite 
number of degrees of freedom, one obtains as starting point a problem of classical 
mechanics. In this way, the basic difficulty lying in the structure of the Lagrange 
equations of motion (being of the first order) can be solved on a purely classical level, 
independently of other problems of essentially quantum origin (as for instance that 
of connection between spin and statistics). 


II. The classical Hamiltonian formalism 


We assume a classical dynamical system of 2f degrees of freedom, with the gen- 
eralized coordinates divided into two groups 4), 9; (7 = 1... f). Let the Lagrange 
function be of the form 

L= An dj 9 — А49 — Kin 9 4 (1) 
(Einstein's rule of summation over repeated indices is assumed; j, k,.... run always 


from 1 to f). The form of L is not so general as given by Schwinger (1953a), but gen- 
eral enough to allow a subsequent generalization to the known cases of fields obey- 
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ing first order equations, e. g. the Dirac or the Kemmer equations. The Lagrange 
equations resulting from (1) are 


2 Asi d; - Ka q,— 0 


J 


245, d, + Кыа, = 0, (2) 
According to the usual definition of the generalized momenta, one obtains 
9L н 
DS Д. 4; 4% 
ЖЕЛДІ: ; (3) 
PES 24 х= 27542 


Nevertheless, а transition to the Hamiltonian formalism is not so immediate іп this 
case. Ав the Lagrange equations (2) are already of the first order, the introduction of 
р} р; as independent variables besides 4» 9; must make half of them superfluous. 
Indeed, the relations (3) have the form of additional conditions between the p's and 
45. А way of avoiding this difficulty (being a classical analogue of Schwinger's quan- 
tum-variational considerations) were to seek for the true, independent pairs of canon- 
ical variables (whose number must be f and not 2/) satisfying equations of the Hamil- 
tonian form. But one.can also investigate whether (or under what conditions) all 4f 
variables 4» d; р} р; сап be formally used as canonical variables, in spite of their 
mutual dependence given by (3). This possibility is suggested by the existing meth- 
ods of treating similar cases (it suffices to mention the well-known problem of the 
canonical formalism for the electron field). Following this second way, one must 
find a Hamiltonian formulation for the variables 4» 4» Рр р; which is compatible 


with (3). The best way to ensure this compatibility is to treat the relations (3) together 
with their solutions 


qr = Cy pi 
dk er АР» - (82 
where | 
Aj Cy, == А, Cu qp д : (3) 


+ 


as а canonical transformation independent of time. This is evidently possible f the 
matrix A is not singular (which we assume in the following). As a necessary condi- 
tion for the correctness of this treatment, the invariance of the Hamilton equations 
and the Hamilton function under transformation (3) and (37) must be proved." 

Adopting the usual definition of the Hamilton function Н == р; d; + ра = 15 


one obtains, using (3), 
H = Kj; % (4) 
eB E PLC ee a с ы re E e ss | 


1 The general postulate of invariance of H becomes the postulate of invariance of its form as, after 
such a transformation, H remains a function of the same 4f variables. 
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(the terms with ў; 4; cancel each other). This is not a unique way of expressing H. 
Owing to (3) and (37) it can be brought into many other forms; but only one of them 
fulfils the postulate of invariance under these transformations. Denoting this *symmet- 


rized by means of (3), (3’)” form the of Hamilton function by H, we have 
Л ЖЕ еу a v a, 
Н = T GC, 4 + Kj, Gi Ve P1 — Kj Cy G pj — Қ Ca Cu Pm pr) = Kir G qr (4) 


where q;, q; are the symmetrized expressions 


а TOR ТЕ 
ia (qj + Сар») 


ЕЛЕЕ (5) 
dj = (aj — Сыр». 
In the same way the р;з can be symmetrized 
ESEP AN A 4 
Bj = 5 (Pi — Au qs) 
zur s (5”) 
Pi = Gi + Ai 98); 
(these formulae will be used later). ң 
It is easy to verify that the equations 
т dH . ӘН 
EE Жар 
ан он (6) 
= gp — зр 


аге invariant under (3), (3’) considered as а change of variables. One must only take 
into account the relations 


9 9 9 9 
ее р р Е 
д; | Әр, Әру i 9q, 
9 д 9 д 
co md Жол kj 777» (2 
94; 9p, Әр; Әр, 


valid for an arbitrary differentiable function of the 4/5 and p's, symmetrized in the 
above sense. The explicit form of the equations (6) with H given by (4) is 


a! 1 ” 
Pr = — = Кы 
a le 2 
kou qi 
2 (8) 
E 1 - 
ф = — 5 Сы Ки 


.” 1 ” 
qk = > “ik Kj qi 
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(8) together with (3) are equivalent to the Lagrange equations (2). Then one can consid- 
er (8) as Hamilton equations, and 4» 9, р} D; as canonically conjugate variables. 
It is to be emphasized that all variables should now be treated formally as independ- 
ent, their mutual dependence expressed by (3) or (3) being already taken into 
account by the symmetrization procedure, which changes the róle of these relations 
into that of a canonical transformation. To show the consistency of the above version 
of the Hamiltonian formalism, a discussion of the Poisson brackets is required. When 
F and G were two arbitrary functions of 4f canonical variables 4» dj р» pj, the 
usual definition should give for their Poisson bracket 


IFIG РӘС | 98F8G OFC 
[F, С] к = + ES 2 M T 7 E cae ” nio (9) 

qj p; 94 Әр) Әр; 94; Әр; 94; 
(Тһе index (с) is introduced to avoid confusion of the classical Poisson brackets with 
the commutator FG — GF, introduced in the next chapter). However, in the presence of 


the supplementary relations (3), the unique sense of the symbols F, G, and their par- 


tial derivatives with respect to the variables q, p can be guaranteed only by the same 
symmetrization rule as used before (for H), with the simultaneous treatment of all 
variables as if mutually independent. At the same time, this is a sufficient condition 
for invariance of the Poisson brackets [F, G] under the transformation (3), (3^), as can 
be verified by means of (7). This shows that, within the above formalism, the notion 
of the Poisson brackets can be used consistently for arbitrary functions symmetrized 
by means of supplementary relations between these variables treated formally as 
canonical variables. 

In particular, on replacing successively F and G by the variables йд» %» Б» B; 
(the symmetrized form (5), (5) must be used here as the canonical variables play 
now the róle of functions introduced into the Poisson brackets), one obtains 


e UNE 1 
[gr Piko = [Ges Ви] = 3 бы 


ПІ 0 1 

[qa gi Xo em 2 Сы (10) 
ae 1 

[Pa Pilo = — 5 4n 


[дь dılo = Ee dro = [Be Plo = [Pe Piko = Ve Pele = Eie Pilko = 0. 
The appearence here of the factor 1/; is characteristic for the case of first order equa- 
tions of motion. Owing to this factor the Hamilton equations (8) result from the gen- 
eral operator equation 

SEALE (11) 
as may be verified by equating . F to jj. qj» ру, апа р; successively and using (10). 


The identities g= = = y; g = = 4» р; == = р, pj шер; allow us to eliminate the symmet- 


rized expressions after performing all the differentiations with respect to these varia 


` bles. 
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ІП. The canonical commutation relations 


According to the well-known correspondence rule of canonical quantization, one 
has to replace the Poisson brackets by commutators or anticommutators of variables, 
treated now as Hermitian operators: 


> Р Я [F, С] = FG - СЕ for boson statistics (12^) 
th [F, С] > Ө (Ё Су е Ж (ЕС un GF), Ө = +1 for fermion statistics (12”) 


(е) 


Applying this procedure to the equations (10), we get the canonical commutation 
relations, at first for the symmetrized operators 4 and p; but making use of (5’), (5”) 
we can also obtain analogous relations for the g's and p's, namely 


[qe Pil = [дь pi] = 3 бы 


[д qi] E AS: Сы (13)) 
[рь pil id zx An 
[ge qi] = [дь 4] = [Px Pil = [рь pi] = [дь рі] = [дъь pr] = 0 
for the case of boson statistics, and 
(a pi} = – Ө а ды 
fan pi) = 0 5-6 
(13^) 


TA а) -- ر‎ 5 Сы 


е, і 
{Pr pis = —@ 2 Ap 


{чь qı} = {дь qi) = (p pi} = {рь pr) = (ae Pi} = (ae p) = 0 
for that of fermion statistics, with Ө =+ 1 (Units ћ = c = 1 are used). Thus the 
necessity of symmetrization becomes superfluous after the passage from the ,clas- 
sical“ to the „quantum Poisson brackets“ since this last does not involve differentia- 
tions. The commutation relations (13^) or (13”) are in fact not independent, owing to 
(3), (3°), but, as is easily seen they, are mutually consistent. Their correctness can be 
verified by equivalence of the operator equations 


with the equations of motion (in the Hamiltonian or in the Lagrangian form). The 
proof, using (13”) or (13"), is similar to that given beforehand for the classical case, 
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with the only difference that now the Hamilton function does not need to be symmet- 
rized by means of (3), (372. 

Thus we have shown the possibility of treating the variables Dj р; as canonically 
conjugate to 4» 4; also in the quantum formalism of the system described by (1). 
The correspondence method of obtaining the canonical commutation relations re- 
mains generally valid, as the characteristic factor !/; appears alread- in the classical 
theory. 


IV. Dependence of the canonical commutation relations on the form the Lagrange 
function 


It is well known that (137 ог (13”) are not the only possible forms of the canon- 
ical commutation relations for a system satisfying the equations of motion (2). 
These equations are invariant under adding to the Lagrange function a time deriv- 
ative term: 


m d ГА ” 

L=L+ ZR W (qj, qj). (15) 
ӘЛЕ b, YE 
9dj' P 9%” 


In consequence, the canonical commutation relations cannot remain unchanged. It 


but such a term will change the explicit expression for p; — 


seems at the first sight that an infinite variety of canonical commutation relations, 
attached to different forms of the Lagrange function, will exist for the same equa- 
tions of motion. However, to be right in calling рур; the canonically conjugate 
momenta, one must investigate whether in all cases their dependence on q}, 4, has 
a character of a canonical transformation. 

All Lagrange functions leading to (2) can be expressed in the form 


=: а 
e (А» 4; дь x Ак 4; qr = Kj, dj qr RS LACE q е (16) 
where E 
W (4,4) = Wd 4» (16') 
or shortly т 
= a (aj, 4; 4, — bin d; Ik — Kir d; 4)» (17) 
with 
Vik = 4; par Wir 
bin = Ay — Ға (17) 


(for generality an arbitrary multiplicative real constant а >> 0 is also introduced). 
From (17) one obtains 
pj = аа 4% (18) 


? The value +1 of © remains undetermined (but the same in all formulae (13")); it depends on the 
order of the operators di: q; in H. 
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and z 

H = pj ij + pj qj = L= a Kag de (19 
A necessary condition that (18) be a canonical. transformation is the condition of 
solubility for 4,4, which is evidently not fulfilled if accidentally 4j, = + Wy. 


These two cases must be treated separately. 
(i) Assuming 4 = W;,, one obtains (denoting this case by the superscript (--)) 


LÐ = а (245, d; Ue — Къ 9) | (20) 
, ILH 0 
PRICE = Ea 21 
J 94; | ( ) 
” oL 
KG) = Dig "y 
Dj 8j; jk qr- (21 ) 


The elimination of one half of the variables is here immediate, as p? vanish identi- 
cally while q; can be expressed by p; «(+ owing to (21”) (the matrix А is not singular); 
the remaining variables qj, p; «^? are already independent and therefore canonically 
conjugate in the usual sense. The Hamilton function expressed by them reads 


1 m m 
Н = PES Capi дь (22) 


and, as can be easily verified, the Hamilton equations 
OH +) 


HE ` 
b; a Q3) 
" OH) 
qi 8p; P 


are equivalent to the Lagrange equations. The canonical commutation relations have 
in this case the traditional form d the factor 1|,) 


[4 pr | = ióy 
(4, 41 = Ір,С9, pj ?] = 0, (24^) 


(disp?) = Ө іды 


| {Чь а} = {Pe Pi} = 0. (24) 
(ii) Quite analogous results сап be obtained for 4, = — Wi: 

I? = — a {2 Aud; de + Kir 9) (25) 

да ILO) > 

qr x 
j = — 2a Áy , 
Di ad kj qk >, (26°) 
Ж ILO) 

pj = ———=0 (26”) 
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and after eliminating que 


1 RT 
HO = — Қа Cn ар” (27) 
(С) aH c 
Di 7 (287 
99; 
v oH 
ques ЕТЕУ (28”) 
Мер °] = id 
[дь ql = [007 рК} = 0, (29) 
ог 
{de PK} = i O ды 
{дь 41} = {рК p} = 0. | (29”) 
For А, F+ W;, equations (18) can be solved for 4,, 9; 
Jum L е ” 
Ч x kl PI (18) 
; 1 =; 
qe = — = Фар» 
where 


ад, сы = ар С = б) 

bj, dy — REO б; dy, = f др. (18”) 
It is also possible to _symmetrize the Hamilton function (and any other function of 
the variables q;, q;, p; ру) using (18) and (18”) treated as a change of variables. De- 


noting this new symmetrization (more general than that defined in Chapter П) by 
the same symbol ~, we have 


= , n 1 1 | о 
Н = 2 Kad goes L Kj сл Tk PI — — Kir dis qj Pi — ai Bik Git dns Pui 


4 
ТЕ ee (19") 
= 4 jk dj ЧЕ: 
However, introducing the transformation (18), (187) into the equations 
REESE enr 
Pii 9g x 0I gps. | (30) 
n^ e H of? 9 Я 
j pj me. Aur E 
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and using the relations similar to (7) 


8 OSG 1.2 
— = QAM. == =й, 
Эа; Әрь Әр; а 9а, 
ә 2, "EE. (31) 
— = a б, Tu ау 8 АО 
29; Op, Өр @ 94, 


holding for all differentiable functions symmetrized according to the new rule (given 
for H in (19) ), one sees that (30) remains unchanged if and only if 


а ан = ар; dy = д, 
5, Сы — б; Ce ді Р (32) 


Іп order to satisfy (32) simultaneously with (18"), one must put 


(33) 
Wa = 0, (33) 
on account of (17^). This last equation shows that only for a Lagrangian which is 


symmetrical with respect to integrations by parts (differring from (1) only by an 
immaterial constant a) 


LO = a (Aid; 4 — 44 4, — Қаға) (34) 
it is allowed to treat 
, ILO n ILO 
| , pj = چ‎ (35) 
9d; 94; 


ав momenta canonically conjugate to dj, 9. This form of the Lagrange function 
(besides the two special cases L‘*?, L) becomes then pointed out among all L of the 
form (16) giving the equation of motion (2). These three cases can be characterized by 
the fact that for them (and only for them) the Lagrange function is equal to zero 
when the equations of motion are fulfilled. 

For the sake of completness let us write down the canonical commutation rela- 
tions for 1/9, which when a 4 1 are somewhat more general than (13^, (13”): 


[de pi | = [dk pi 9] == 5 бы 
"MEAT, L 
9 di] = — 5 Сы (36) 


«0 "(0 ia 
[pe Pi ] = — > Фа 
with all other pairs commuting, or 


i 


fap, p) = = 6 2 


ды 
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ха І 
{qr РІ S =Q D ды 
wr, i 
{ 4 qı ) = 0 о Сы (36”) 


” га 
р.р) - өзе 


with all other pairs anticommuting. 


V. A generalization of the results for the case of field theory 


The above considerations applied to a model with a finite number of degrees of 
freedom. But the results can be immediately extended to the known cases of fields 
described by first order equations, i. e. to the Dirac and the Kemmer fields. Tt suffices 


to exploit the analogies 
Grp dv 
TURAE TIR Иб 
ол CEDE — 
Dj oy Р; Bu (37) 
i 2 
A = ғ Vu [o 9» Й 
С + — 21у, (ог — 2if,) 
К >x, 


y, B, denoting Dirac and Duffin-Kemmer matrices respectively, and x being the 
mass of the fermion or boson particles (in units ћ = c = 1). It is easy to see that, 
using (37), one obtains the well-known expressions for the Lagrange functions (sym- 
metrical with respect to integration by parts) 


1 E БЕ 767 
дов ду eB yoy) ea. (38°) 
or | 
1 ==. -- -- 90/1 
£0 = af 3859. - o9 fey) — ev. (38) 


leading to the Dirac and Kemmer field equations respectively. Analogously, the 
commutation relations (36^, (36") allow a straightforward generalization 


(р по} = Ө 5 д, Ó (x — x’) 


ѓ (14—05 beo (e = x) 
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(ро, Vo) = (9 E (V4)ec à (x E х) (39) 


{ль Tic) = O F (уде 6 – =) 


{Wop Pod = (ро Wo} = (п Ta} = (ло по} = (onm) = (ұр по} = 0, 


or 


[vo лә] = [os] = = ды Û (% — = 


E. vd] = — — (80 5 (x — =) 


cS (39") 
[n 32] = — 2 Bader Öl — x) 


[vi val = [рь vel = [p ло] = [ль ла] = [yo ла] = [ve ло] = 0 

for the two fields respectively (with t’ = 2). In (39") the index I denotes the dynam- 
ically independent components of the field variables y, y that remain after elimi- 
nation of the dependent ones (which appear in consequence of the constraints for the 
Kemmer field); this procedure excludes at the same time the eigenvalue f, — O (for 
details and references see Hanus (1954)). 

The form of the relations (39) suggests the use of a normalization convention 
a = 2 supplying all the formulae with the same factor !/, (although this factor is then 
of a different origin in the first two and in the following two formulae (39). 


КРАТКОЕ СОДЕРЖАНИЕ 


В. Ганус, О формализме Гамильтона и канонических соотношениях коммута- 
ции для случая уравнений Лапранжа перво1о порядка 


В настоящем сообщении для случая уравнений движения первого порядка, 
рассматривается проблема формализма Гамильтона, на модели динамической 
системы с конечным числом степеней свободы. Показано, что можно ввести пере- 
менные сопряженные канонически с обобщенными координатами согласно с их 


обычным определением р; = Tks нестмотря на TO, что они зависят OT переменных 


aL 
94, 
если эти зависимости рассматривать как каноническую трансформацию незави- 
сяшую от времени. Каноническое квантование при помоши принципа соответ- 
CTBeHHOCTM (использываюшего свойства скобок Пауссона) можно провести фор- 
мально тем же самым способом, как и в обіцем случае; получаются правильные 
канонические правила перестановки содержащие фактор 1/» Возможность pac- 
сматривания соотношений между переменными Q и р как каноничной трансфор- 
мации ограничивается к одной форме функции Лагранжа из всех бесконечно 
многих ведущих к тем же самым уравнениям движения (т. е. отличающихся 
между собою членами являющимися производными по времени любой функции 
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переменных 9); формой этой является форма L симетрическая в отношении 
к интегрированию по частям. Следовательно эту форму (связанную с правилами 
перестановки содержащими фактор 1/2) надо считать особенной; кроме того отли- 
чаются от иных также две другие специальные формы, разрешающие на не- 
медленное устранение половины переменных Q. Остальные не зависящие KaHO- 
нические переменные выполняют тогда традиционные правила перестановки (без 
фактора !/). Результаты разсуждений были обобщены на известные случаи полей 
подчиняющихся уравнениям поля первого порядка (поле Дирака и поле Кеммера). 
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SIMPLE MODEL OF THE ‘Li NUCLEUS 
AND THE NEUTRON INDUCED REACTIONS ON *Li.*) 
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A simple model of the 9Li nucleus consisting of an a particle plus a deuteron is 
discussed. This model is applied to the calculation of the differential cross-section for the 
SLi (n, 4)“Не and Li (n, d) He reactions. For the first reaction the pickup mechanism is 
assumed, for the second one the mechanism of a „knock-out“ of the deuteron as a whole 
from the *Li nucleus. The results are compared with experiment. Fairly good agreement 
is obtained especially for the 9Li (n, t)4He reaction for small angles. 


8 1. Introduction 


On interpreting their experimental results concerning the Li (y, 4)“ He reaction 
Jensen and Gis (1953) suggest a simple two-body model of the Li nucleus consisting 
of an a particle and a deuteron. This assumption was supported by Vashakidze and 
Chilashvili's (1954) theoretical results. 

In the present paper the Li (п, t)* He and 614 (n, d) He reactions are investigated 
theoretically from the point of view of the „a + d^ model of the Li nucleus. A strip- 
ping-like mechanism is assumed for both reactions. The Li (n, t) “Не reaction is con- 
sidered as a pickup process of the deuteron as a whole and the 9Li (n, d) He reaction 
as a ,knock out“ of the deuteron as a whole by the incident neutron from the 614 
nucleus. The results are in fairly good agreement (especially those for the 9Li (n, t) 
4He reaction) with the experimental data (Frye 1954, Weddel and Roberts 1954). 

It is well known (Gerjuoy 1953) that all the widely used methods in the up to 
date stripping theory may be reduced to the Born approximation. To obtain expres- 
sion for the differential cross-section in the form given by Butler (1951) or Bhatia et al. 
(1952) some additional simplifying assumptions are, however, necessary. In the pres- 
ent paper no additional assumptions are introduced because of the knowledge of the 
wave function for the model involved in our treatment. This knowledge gives the 
possibility- of performing the exact Born approximation similarly to what was done 
in the theory of the ?Be (p, d) 8Be reaction (Dabrowski and. Sawicki 1954, 1955). 


*) This is a further development of a paper by the same authors (1955b). 
(323) 
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In $2 the assumed model of the SLi nucleus is discussed. In $ З the theory of 
the Li (n, t) “Не reaction is given and compared with experimental data.,In $ 4 an 
attempt of the explanation of the Li (п, d) He reaction and its comparison with exper- 
iment is given. In the Appendix the formula for the Born approximation cross-section 
for the Li (n, t)4He reaction is derived. This derivation, although not identical with 
the standard Green’s function technique, contains nothing essentially new. For the 
Li (n, d) °Не reaction the derivation is similar. 


82. Discussion of the „а + d“ Model of the “11 Nucleus 


The model of the Li nucleus consisting of an a particle plus a deuteron is support- 
ed by the following reasons: 

1) The (y, d) threshold energy is very low (1.477 MeV), the average binding 
energy per nucleon being — 5.3 MeV, the threshold energies for (y, n) and (y, p) 
being respectively ~ 5.8; 4.7 MeV and the energy for the disintegration Li -+ y = 
= 3?H (= ЗН + ЗНе) ~ 25 MeV (~ 16 MeV). 

2) The total cross-section for the (y, d) reaction is extraordinarily small (<3.5 - 

10-29 cm? for 2.62 MeV y-energy—Jensen, Giss 1953). 

3) The spin I — 1 equals to the spin of the deuteron, while the a particle is spin- 
less. 

4) The magnetic moment u = the magnetic moment of the deuteron. 3) and 
4) give us the reason for the S-state assumption for the ,,a + 4“ system). 

5) The present model was satisfactory in the explanation of Glenn's (1952) 
and Jensen & Gis's (1953) results for the (y, d) reaction (Vashakidze and Chilashvili 
1954). 

6) It is also consistent with the Lauritsen, Huus, and Nilsson (1953) results con- 
cerning the a- d scattering. 

The assumption of a stripping-like mechanism for the Li (n, t) “Не and Li 
(n, d) *He reactions is the consequence of the assumed model of the Li nucleus. The 
immediate evidence of this mechanism is given by Frye's (1954) experimental results 
for 14 MeV incident neutrons. The total cross-sections for both the reactions being 
comparatively large are respectively (26 + 4) mb,(89 +5) mb (the total cross section 
for the SLi (n, p)8He reaction (is 6 + 2) mb) and tritons (deuterons) are scattered 
mainly in the forward direction. 


63. The SLi (п, t)*He Reaction 
(A) Calculation of the Differential Cross-Section 


In the present calculation the Coulomb effect is neglected. This assumption is 
justified for larger neutron-energies, e.g. for 14 MeV neutrons applied by Frye (1954). 
The interaction V,, between the incident neutron and the a particle is also not taken 

5; М 


1 The binding energies were calculated from the data given by Li et al. (1951). 
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into consideration. In the case of the stripping reaction Horowitz and Messiah (1953) 
have shown this interaction to give some decrease of the absolute values of the cross 
section without any significant change of the angular distribution. 

Under these assumption the Born approximation result for the differential 
cross-section in the c.m. system is (see Appendix): 


do 1 
10 — Gam М My 573 za (1) 
ошаш, 
where 
I= f dî, d£ dù d$ do, do, do, er, X}, XG V4, Po (2) 


where: M; = (12/7) M, М; = (6/7) M are respectively the triton and the incident 
neutron reduced mass (M being the nucleon mass); hk,, hk, are respectively 
the triton and indicent neutron momenta in the c.m. system; ро, U, pg are respectively 
the neutron, triton and Li nucleus magnetic numbers; X, = Х.(2) is the internal 
wave-function of the a particle (& being the internal coordinate of the | a particle) ; 
Хи, 1 
(ғ, -- т„)/2, и = \3/2 (7, = 7p) are the internal coordinates of the triton; Fy, r,, 


is the internal function of the triton; 075 are the spin variables; = To 


r, are the vectors from the center of mass of the Li nucleus to the incident neutron 
and to the neutron and proton of the deuteron in the Li nucleus (which is picked 
up by the incident neutron) respectively; 7, is the radius vector from the a particle to 
the triton; Vja is the interaction between the deuteron and the a particle within the 
Li nucleus, and 


Чо = AT (00) ейт I (3) 


is the wave function of the incident neutron and the 611 nucleus in their initial states; 
X,» Xp, are the internal wave functions of the neutron and the Li nucleus respectively. 
In agreement with the assumed „а + d^ model of the Li nucleus we put 


Dn (m, o.) E ie u) Ф (raa) Xa (2), (4) 


ү2 
Il 


1 : 
а шш) Xun (0,) ° Xw, (0p) is the deuteron 


4 


where S, (Op 6,) = 2 2 Hull 
“лир 
spin function (Хи,, Хи, are respectively the neutron and proton spin functions and 
LI 1 

Ga Ps CHE ET 

deuteron function of the 1S state (the 1D state admixture is neglected) and Ф (т) is 
the internal function of the „а + d^ system, y,, being the radius vector from the а 
particle to the center of mass of the deuteron. 


2 
n; is a Clebsch-Cordan coefficient); У m үз“ 1 u) is the internal 


25 | 
For YW (ғ u) wë put the unperturbed internal deuteron function. Thus we 
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neglect the distortion of the deuteron in the Li nucleus. This fact seems to be the least 
satisfactory assumption of this calculation, but this assumption may be supported by 
the following consideration: The incident neutron picks up the deuteron as a whole 
from the Li nucleus only when at the moment of the collision this deuteron is undis- 
sociated. However, while the neutron and the proton of the deuteron are within the 
range of their mutual interaction, the perturbing interaction from the a particle dis- 
turbs the п-р relative motion comparatively weaker. So we may hope that the distor- 
tion of the deuteron in the 614 nucleus is of no essential importance for the 9Li (n, t)*He 
reaction. 


For Vj, we assume Гм = Vaa (Taa). 
For X,, we put 
Ха, = 5и, (Fp 6, 0o) В, (s, и), (5) 
where К, (s, u) is the S ground-state function of the triton which is symmetrical in the 
spatial coordinates of three nucleons, and 


1oo) xe 90 (6) 


ЕГ 
Su, = Хи, (0p) У (1 5 — Шо Шо 


is the triton 2S —state spin-function antisymmetrized in the spin coordinates о” and o, 
of the neutrons. 

Since all the functions in (1) are known, the integral in (1) may be axactly calcu- 
lated. Making use of the equations 


ry = (2/3) r, — s and r, = г, — (1/3)"s, we obtain 


1 ` 
58», ИВ = S| Tel Ik |, (7) 
Бойы 
where 
1 : ena 
S = 2:3 У | da, do, da, Sut (бр On 00) Sy, (Gn Fp) Хи, (оо) | = At (8) 
Moe e 
I, = | dra e V, o o (Pad) (9) 
+ J+ „iks 2 
1, = | е В; (s, и) У 5“) > (10) 
where | 
Ё 2 > ruga 4 2 ۾‎ 
RET ES p tog ko + kr — y ko Е, cos 0 
e DRE ES 1 2 11 
= —№ k? = ko + > kr №0050, 9 


0 being the c.m.s. scattering angle. 
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On assuming for V}, а square well of radius Ry and depth У, we have 


7-Sln X Tad 
— forira T 


Ф == Tad (12) 
Wi e Тоа е-аа- в) for Tag < Ro, 


Веера = M З Mahe , = 2 Mi - Th, А = - ar (1+ aR) 


Center of mass diff cross-sect. 
(n millibarns / steradian do/da 


e о 


40° 80 120 
Center of mass angle @ 


Fig. 1. Angular distribution of tritons from “Li (n, 1)“Не for 14 MeV neutrons. a denotes the exp. 
curve; 6,c— the theor. curves (relative values). for Ко = 4۰10713 cm and 6۰10713 cm resp. 


do/do 


Center of mass diff. cross - sect. 


іп barns/steradian 


0 40° 80° 120° 160° 
Center of mass angle @ 
Fig. 2. Angular distribution of tritons from 614 (n, t) “Не. a, d denote resp. the exp. curves for 1.5 and 2 
MeV neutrons;, b, c — the theor. curves (relative values) for 1.5 MeV neutrons resp. for Ry = 4: 
- 10-13 сп and 6: 10713 cm; е, f — the theor. curves (relative values) for 2 MeV neutrons resp. for Ro = 
= 4۰10713 cm and 6: 10713 cm. 


(= = 1.477 MeV being ді binding energy of the deuteron in the SLi nucleus.) 
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А ee 


For I, we find 


Тк = 4 


х-К x+ К 
For R, we assume two forms: (1) R, = (2y)3/r le 290049 (a gaussian function), 
where у = 2,5 - 1012 стт! (Newns 1952); (2) R, = № e-V2« Vus [(u? + s2)", where 
N,is the normalization factor and 1/a = 1,318 - 10713 cm, n = 1/4 (Irving 1951). 
For the deuteron function W we assume two forms: 
а) V/(o) = (30/4л) “е” “440, where ag=VMe,|h?, eq = binding energy of the deuteron 
(the exponential function with the normalization of Smorodinski (1948)); b) ф(о) = 
= М (её е 1%) 0, where Ny = (1/22) |8Г(8 + Г)]/(р — Г)%Г = 1,234: 10cm, 
В = IT! (the Hulthen function with the parameters given by Kruse et al. (1953)). 
The integral Г, is discussed in the three following cases: (la) and (15) (К, being 
the gaussian function). In these cases 


sin («—K) Ro ee (13) 


co 


Т == 8 V2x | [* u? P (% u) e-2y joe (14) 
0 


The angular distribution is thus independent of the choice of the form of the deuteron 
function. Only the absolute value of the cross-section depends on the deuteron function. 
The integral in (14) may be immediately transformed to the probability integral for 
both kinds of W ((a) and (5)). (2а) (R, being the Irving's function, W — the gaussian 
function). In this case we have 


1, = > laf NÉT (4—2n)/ (2a2)?-7. С (5) ; (15) 
where 
1 
С (р) = 2 f dx sin [(4—2n) tg-! B/A)] | (В? + A2)2-", (15”) 


0 

where 4--1--(2/3)% (а, /а)х, В = (kjay2) (1 — x2)^. Тһе angular dependence 
of Г, is connected with the factor С (p), the graph of which was given by Newns (1952). 
With the help of this graph we find that the angular distribution in case (2a) differs 
perceivably from that in the cases (1a).and (1b) only for larger neutron energies and 
for large scattering angles (some increase of do | dQ for large angles is obtained). 

For case (2b), i.e. for Irving's function К, and Hulthen's function ¥, the depen- 
dance of J, on k is similar, but somewhat more complicated than that given by Eqs 
(15). So we omit it in the comparison of our theoretical results with the experimental 
data. | 


(B) Comparison with Experiment 
For Rg we use the value Ry = 4 - 1071? cm suggested by Lauritsen et al. (1953) 


(the corresponding value of V is V = 5.772 MeV) and the value Ry = 6 - 10^ !3cm 
which lies in the region of values of the ,,channel radius“ for the elastic a — d scat- 
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tering discussed by Lauritsen et al. (1953) (the corresponding value of psp es 
= 3.679 MeV). 

The resulting angular distribution for the cases (1а) and — (1b) compared with 
Frye’s (1954) results for an incident neutron energy E = 14 MeV is given in Fig. 1. 
Although the theoretical values of 4 /dQ for large angles are far too low, as usually 
in the stripping calculations (due mostly to the compound-nucleus background), we 
observe a general agreement of the shape of the curves. At large angles we obtain 
a second maximum, contrary to the case of small energies (Fig. 2). This maximum 
for Ry = 4- 1071? cm is scarcely visible in the scale of Fig. 1. 


Table I 
Theoretical values of (40/0) 0 = біп mb/ster for the Li (n,t) “Не reaction. 


Irving’s 


58,1 54.4 11,6 
23,0 24.5 1.3 


Table II 


` Theoretical and experimental values of (45/40) 0 = 0 for the Sri (n, “Не and 614 (n, d)|He reactions 
for 14 MeV neutrons. (For (n,t) the best values of $3. are given). 


Theoretical values of (40/40) іп mb/ster. Exp. values of 
(40148) — о іп mb/ 
Ry = 4.10 cm | Ry = 6.10 cm |№ = 7,5.10 ??cm| ster (Frye 1954). 


38 
7.5 


Li (n, d)$He 78.8 
Li (n, t))He 


However, for smaller neutron energies, when the Born approx«mation .seems 
to become more doubtful and a distinct Coulomb effect appears, it se ems to be also 
interesting to compare our results with the results of Weddel and Roberts (1954). 
From Fig. 2 we see that a better agreement with Weddel and Roberts' results for E — 
= 1.5 and 2 MeV was obtained for Ry = 4: 10713 cm. Since in stripping the Coulomb 
effect gives some flattening of the angular distribution curves (Yoccoz 1954, Sawicki 
1955 ab), we hope that by taking the Coulomb effect into account the agreement of 
t he theoretical and experimental curves might be improved. 

The angular distribution for case (2a) is nearly the same as for the cases (1a) 
and (1b). Only for 14 MeV neutrons case (2a) gives a slight increase of do/dQ for larger - 
angles, as already mentioned in the foregoing section. 

The absolute values of the calculated values of (do/dQ),.., are given in Tablel. 
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The corresponding experimental values аге: ~ 46 mb/ster for E = 1.5 MeV, ~ 27 mb 
ster for E = 2 MeV, ~ 7.5 mb,ster for E = 14 MeV. From Table I we see that for 
the gaussian function R, the use of Hulthén’s function Y increases (do[d£2), —о 
over 2 times the respective results with the exponential function W. Consequently for 
R, = 4: 1073? cm no better agreement of the absolute values of the cross-sections 
with the experimental values is expected in case (2b) but for Rọ = 6: 10713 cm the 
agreement in this case can be fairly good, especially for E — 1.5 MeV and also for 
E — 14 MeV. 


84. Тһе SLi (n, d) “Не Reaction? 
(A) Calculation of the Differential Cross-Section 


In this section we give an attempt of the explanation of the 9Li (n, d) He reaction. 
Since the calculations for this reaction are somewhat similar to those of 83 we give 
less details. 

Under the assumptions of $ 3 (here the neglect of V,a is equivalent to the neglect 
of Vja because of the equation similar to (7A)), we obtain for the Born approximation 
cross-section in the c.m. system: 


do 1 EET ; 
Lp ДР ажаты 16 
40 бақа» Ms Mo org 3, UIS (о) 
y 
where 
I — (d£dr,, do dr, doy do, da, X; Ху еба "4 У, Po, (17) 


where M; = (10/7) M is the He — deuteron system reduced mass; ћ Ё is the deu- 
teron momentum in the c.m.s., Хи is the ?S-state deuteron function, wz being the 
magnetic number of the deuteron; 


tn = У [Поют у (Bopa) ta (б) Вод Q9 
us my 
is the P5/,-state internal function of the Не nucleus (treated as the „a + n,“ system) 
| being the magnetic number of He; " is the radius vector from the a particle to 
the incident neutron; ғы is the radius vector from the ӘНе nucleus to the deuteron, 
0 E V3u; V a is the interaction between the incident neutron and the deuteron in 
the Li nucleus. All the remaining symbols are explained in § 3. 
We assume the simplest form of V,,, namely 


ӛламһ2 „> 
; Vea = a) (19) 


— the Fermi „contact pseudopotential“ with а, = V(2/3) а? + (1/3]a2 = 5.08. 


® The results of this section were already published in Bull. Ac. Sci. (Dabrowski, Sawicki 1955а). 


pa 
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: 1073? cm — the n—d „scattering length“ given by Christian and Gammel (1953), 
This interaction is consistent with the assumed model in which the deuteron is treated 
as an individual particle. This corresponds to the simple „knock-out“ mechanism of 
the whole deuteron by the incident neutron, the latter being captured by the a particle, 
similarly to the mechanism of the (n, р) reaction proposed by Austern, Butler and McMa- 
nus (1953). Evidently, it is a crude assumption since the deuteron is a fairly loose struc- 
ture. In more exact calculations Z,, should be used in the form of a sum of interactions 
between the incident neutron and each of the nucleons of the deuteron separately. 
In this way besides the „knock-out“ mechanism we would have the pickup of the proton 
by the incident neutron. However, the contribution of the latter mechanism to the 
angular distribution is decreased also by leaving out of account the contribution 
coming from the incident neutron — a particle interaction Рр, By allowing for the 
latter interaction, we would arrive at a pickup mechanism of the a particle from the 
SLi nucleus by the incident neutron which would cause an increase of the number of 
backward scattered deuterons, while the ordinary pickup gives an increase of the 
number of forward scattered deuterons. 
Under our assumption we have after all integrations and summations, 


а 2 3 n2 
IO zs Moz (еі LEE | | | drr? Alan RA | (20) 
where j, is the spherical Bessel io Hon of order 1; 
а? = 204 ee v ны ky cos 0, j (21) 


0 being the scattering к in the c.m.s. 
Тһе «He nucleus radial function R we assume to be 


forr <a, 


жс ыы s 
ТЕРА АН (22) 


Изя for r > а, 
r * 


appropriate to a particle in a P level of zero total energy in a square well of width a 
used by Hughes and Le Couteur (1950). The value а = 2.9 · 10713 cm is used as 
suggested by Adair (1952). 


(B) Comparison with Experiment | 
The resulting angular distribution for 14 MeV neutrons was obtained from (20). 
with the help of graphical integrations. 
——————Á—————SMÁ—— Á—ÓM—Ó————— n— x ——À— тт т” 


? The assumption of zero range interaction for neutron energy of about 10 MeV must be чыч 
as а rough approximation ‘only. 
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The absolute values of (do/dQ),_, for different values of Ry for both reactions: 
SLi(n, d)»He and 9Li(n, 2)“Не compared with the respective experimental values are 
indicated in Table II. 

While for the 8Li(n, t)*He reaction the best angular distribution was obtained 
for Ry ~ 6 · 10 13 cm, in the present case we get a better agreement with experiment 


in whb/ster. 


observed (Frye 1954) 


30 


ao 
ай 


-73 . 
Calculated (Ry 15-10" "Cm)-retative values, 


20 


Center of mass diff cross- sect. 


0 40° 80° 120° 160° 
Center of mass angle @ 


Fig. 3. Differential cross-section for the Li (n, d)5He reaction for 14 MeV neutrons. 


for larger values of Ry. The angular distribution for the 9Li(n, d)*He reaction for Ry — 
= 1.5۰10 18 cm is given in Fig. 3. For smaller values of Ко the maximum of the 
curve becomes broader and somewhat displaced towards larger angles. 


We see that in spite of the simplifying assumptions a general agreement of the 
characters of the theoretical and experimental angular distributions was obtained 


for the Li (n, 4)5Не reaction without taking into account the „ordinary“ pickup 
mechanism. 
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APB BENDIX 


Deduction of the Born-Approximation Cross-Section for the Li (n, t) He Reaction 


The calculation is performed in the c.m. system. The explanation of all the sym- 
bols used is given in $3. All the vectors used are indicated in Fig. 4. 
The Hamiltonian of the system is of the form 


Н- Не(таос„о,) + д i Vat Jes چ‎ H(& T H (su0,0,00) ar а y T7,.(1A) 


Н, is the internal Hamiltonian of the Li nucleus, Н, — the internal Hamiltonian of 


4 


o 


Fig. 4. Situation of the particles involved in the SLi (n, t)*He reaction. Sg, S, , Sg are the centres of mass 
resp. of Ње Li nucleus, the triton and the deuteron; n, denotes the incident neutron. | 


the a particle, H, — of the triton; 7, and Tr; are the relative motion kinetic energy 
operators respectively of the system: incident neutron — Li nucleus and of the 
system: triton — а particle; V;, is the interaction potential between the particles 
(i) and (k). The transformation of H in Eq. (1A) is evident. 
We have to solve the Schrédinger equation of the problem: 
HY = EY (2A) 
in the form ; 
Ч == Е Ч, (ЗА) 
where ¥, is the plane wave of the neutron incident on the Li nucleus (see Eq. (3)) 


and YW, represents the triton departing from the a particle. 
On substituting (3A) into (2A), we obtain in view of (1A) 


(H, а Н, IF 15 НРУ E, — > (Paa E Бо) PECES (Voa c Via) 0 (4A) 
We expand VA, in the following series 


Wy = jas) Gu op o, 00) f. (64) 
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On substituting (5A) into (4A) and making use of the orthogonality of the 7,/8 
and X,'s we obtain for f (7,) an inhomogenous equation, which may be solved with 


help of the outgoing wave Green's function ~ ег. From the form of f (r,) for r,> оо 
we finally obtain the following expression for d |40 


do 1 рі H k, 1 
49 ^ (балар М: М в № 


Вов И; 


[| dr, d£ du ds do, dop do, e^" X, Xa 


{(Vaa "s Voa) H (Voa да Via) Po} (6A) | 


2 


On expressing V,, and Vja in terms of the other operators from (lA) and in view 
of the fact that these operators are all Hermitian and the equations are fulfilled by the 
functions appearing in (6A), one may easily show (Schiff 1949, Gerjuoy 1953) that 


fd, а du ds do, do, doy е-# x, Xa Woa — Va) Ws = 0. (7A) 
After applying the Born approximation we replace in (6A) W by V/, and finally 
obtain 


do 1 ЕСІ 


dO Crh? MeMo „2.3 
> | f dr, dë ай d$ do, doy doy eX, Ka (Vas + Vos) V | (8A) 


Бой 


It is interesting to note that (8A) Vz, may be replaced by V, іп view of Eq. (7A). 
This last result is an immediate consequence of the reciprocity theorem for nuclear 
reactions (see Bransden 1954). 


КРАТКОЕ СОДЕРЖАНИЕ 


М. Домбровский и T. Савицкий, Простая модель arpa) Li и реакции 
индуцированные нейтронами в әтом ядре de 

В работе рассматривается простая модель ядра 6Li состоящего из а-частицы 
и дейтерона. Эта модель применена к численному определению дифференциаль- 
ных сечений реакций 6Li(n,t4He и 6Li(n,d)>He. Для первой реакции принят Mexa- 
низм ,,ріскар”, для второй -- механизм ,,knock-out” дейтерона, как целого, из 
ядра 614. Результаты сопоставлены с экспериментальными данными. Получена 


довольно хорошая согласованность, особенно для реакции 6Li(n,t)4He для не- 
больших углов. 3 
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A generalized field formalism (closely connected with that of Yukawa) incorporating 
a fundamental length / may be regarded as a new type of field quantization: the ,,J-quant- 
ization“, to be distinguished from the canonical ,,h-quantization“. The /-quantization 
appears to be the mass and spin quantization and leads to a rotator (or spinrotator) 
model of elementary particles. The particle is represented by (the surface of) a sphere 
of radius 2М (M being the rest mass). The general solution of the (first or second order) 
field equation includes a whole family of particles with higher spins and masses. A close 
connection with Fierz’s theory of higher spins is discussed. The particles are classified in 
families and super-families according to the type of the field equations (of the first or se- 
cond order) and the type of parity. Four families of bosons and six families of fermions 
follow from general group-theoretical considerations. There exists a homomorphism 
between all single and multi-valued irreducible representations of the group of rotations 
and inversions and the possible types of elementary particles. 

All sorts of particles known at present (light fermions, light and heavy mesons, nuc- 
leons, hyperons) fit into our scheme except the photons (and gravitons) which have to 
be considered separately. Qualitative properties (decay schemes, etc.) are explicable in 
terms of selection rules due to the conservation of angular momentum and parity. Theo- 
retical values of the masses which follow from this formalism agree in some cases remar- 
kably well with experiment. In the remaining few cases (и, л, 2+) the agreement is 
only rough but the discrepancies are probably caused by a ,,fine structure“ of the mass 
levels (self masses) due to interaction. ; 


Introduction 


During the last 20 years an unexpectedly large number of particles in the mass 
interval 200 m, — 3000 m, has been discovered, all of which seem to merit the name 
of elementary particles. 

Under the name „elementary particle“ we understood so far a particle described by 
a field quantity y (x) belonging to an irreducible unitary representation of the Lorentz 
group, and therefore characterized by a mass m, spin ОНЫ, ..., and ‘parity. 

In view of the appearance of more and more new particles with puzzling proper- 
ties the question arose whether we are right in limiting ourselves to group representa- 
tions connected with space-time, or whether we should suitably extend the group 
and the underlying manifold. This last point of view was represented by Pais (1953), 
who regards the problem of isotopic spin as a clue to the problem of new unstable par- 
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ticles. The idea of Pais is very interesting, nevertheless a more orthodox way of 
approach is to investigate first the possibility of explaining the problems connected with 
the new particles within ordinary space-time, before we decide such radical a step as 
an extension of the manifold and the introduction of further dimensions possessing no 
immediate physical meaning. The following reasons favour the more orthodox way: 

(i) The classification of elementary particles with the aid of representations of 
groups closely connected with ordinary space-time (Lorentz group, or its subgroup 
of rotations and reflections in 3-dimensional space, or the closely related unimodular 
group) was already successful in the past. It has provided the fundamental distinction 
between bosons and fermions and the distinction between scalar, vector, pseudoscalar, 
and pseudovector particles. 

(ii) Not all the possibilities of relating particles to irreducible representations 
of these groups have been exploited so far. 

(iii) Until all the possibilities of classification within the ordinary manifold have 
been exploited, an arbitrary extension of the group and manifold leaves too large 
a margin for phantasy. | 

The best known particles (electrons, nucleons, pions etc...) are characterized by 
means of some irreducible representations (single or double valued) of the group 
of rotations and reflections. But it is the group that matters, and not some of its repre- 
sentations, so that it must appear strange why only some and not all irreducible repre- 
sentations of the group of rotations and reflections should be employed in the theory 
of elementary particles. This seems to be an arbitrary discrimination of some of the 
representations from the others, constituting a conspicuous defect of the traditional 
field theory. This defect will be removed if we assume that to every irreducible repre- 
sentation there corresponds a field and if we shall find the lacking types of particles 
among the new unstable ones. Then the situation will be more satisfactory from 
a group-theoretical viewpoint, and we shall possess a clue for understanding the prop- 
erties of elementary particles. 

To fulfil this program a generalized field formalism is needed with a new field 
concept: a field has to describe a whole family of particles characterized by different 
spins and masses. However, the representations of the Lorentz group are characterized 
by discrete spin values but continuous mass values. In order to get discrete mass eigen- 
values the formalism must include an essentially new feature which, most probably, 
springs also from a group-theoretical background. Of course, this may be just the 
group of transformations іп isotopic spin space (and Pais showed the possibility of 
getting a discrete mass spectrum from an eigenvalue problem in isotopic space) but, 
as previously stated, a more orthodox way is to exploit first all the possibilities within 
the common framework of description of physical phenomena. Such a possibility is 
offered by the fact that the Minkowski space enters into physical description twice: 
as ordinary space and as momentum space. The group linking the two spaces is Born's 
group of reciprocal transformations. We shall show a close connection between a mass 
eigenvalue problem and the group of reciprocal transformations. 
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CHAPTER I. А GENERALIZED FIELD THEORY 
l. Some general remarks on the procedure of quantization 


The most essential features distinguishing quantum theory from classical theory 
may be summarized as follows. The manifold of points (ху, X2, x4, %4) = (x, Ee) Aa, 
whose transformations form the Lorentz group, is the Minkowski space. The notion 
of this space will be generalized in a natural way if the quantities x, are regarded not 
as numbers but as operators x, (х аге real operators, x, is an imaginary operator). 
Their counterparts are the infinitesimal displacement operators d, satisfying the 
commutation relations | 

а ор (1) 
From (1) it follows that 4 is an imaginary, d, a real operator. If the time scale is chosen 
so that c = 1, the dimension of dis ст 1. It will be natural and convenient to intro- 
duce, besides c, another fundamental constant J with dimension of length, and inves- 
tigate the dimensionless quantities /71х, Иа. The natural kinematic units are с = 
=з бше р | 

Expressed in natural units the operators x and id are subject to а non-commuta- 
tive algebra. They may be added, multiplied by numbers and one by another. The 
operators x and id span a space which is a generalization of the Hilbert space. We may 
consider canonical transformations replacing the operators x and id by other quanti- 
ties (being functions of the former) preserving the commutation relations (1). In partic- 
ular, we may investigate linear transformations preserving the four-vector character 
of these quantities. An example of such a transformation is 

id, = 5 (eine ld): 


(2) 
iid, = = (і-ік, + ild,) 


Such transformations with two row and column matrices form a group. Later on we 
shall need its subgroup consisting of matrices 


К Же рз (сай es ae 


isomorphic with the group 1,1, —1, —1. It is the group of reciprocal transformations 
(Born transformations). 

From the abstract scheme of operators we may go over to matrices by introducing 
a frame of reference in the generalized Hilbert space. In particular, it may be a frame 
of reference connected with the eigenvectors of the position operators (whose eigen- 
values will be denoted by x). In this representation we have 


0 9 
«x | ار‎ > уду — а), > > = 0800 -x) 0 
qu 
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where ó(x) is the four-dimensional Dirac 6-function. Together with x, and id, 
any operator (being their function) becomes a matrix. 

The quantization of a material point consists in an assignment of the operators x, 
to the physical notion of position in space time and the operators ihd,to the physical 
notion of momentum and energy. In the units c = l = 1, Planck's constant f has 
the dimension of mass. The natural physical units are с = l = h = 1. The assign- 
ment of the operators iAd, to the momentum and energy vector p, is done as follows: 


— ihd is identified with p while the equality between ihd, and p, is regarded as a con- 
dition upon the state vector 

-ihd = p(s) (5) —ihd,|  — pl». (5^) 
(5") is nothing but the Schrödinger equation. In consequence of the weaker connec- 
tions (5") the operator id, does not enter explicitly into the observables and there- 
fore its cannonical conjugate f may be treated simply as a parameter. 

Together with position and momentum, every physical quantity connected with 
the material point becomes an operator, or, in a representation, a matrix. The exten- 
sion of the quantum formalism to the case of generalized coordinates and momenta 
as well as the generalization for systems with more degrees of freedom need not to 
be repeated here. 

We have sketched above the main features of the quantization procedure to 
draw attention to its most essential aspects. In our opinion the most important feature 
of quantization is the introduction of the displacement operators and the fact that every 
physical quantity becomes a function of the position and displacement operators. 
This fact is to be regarded as the general feature of quantization and not the connec- 
tions (5) which are rather of a specific character. 


2. Two types of field quantization 


Let us go over to field theory. The fundamental physical quantities describing 
a field are the field „coordinates“ р(х). Every physical quantity attached to the field 
may be built up of the field quantities and their derivatives as every physical 
quantity attached to the material point may be built up of the coordinates x, and 
their derivatives. The field quantization (second quantization) may be regarded as 
a straightforward generalization of the first quantization for the case of systems with 
infinitely many degrees of freedom. 

However, if by quantization we do not mean necessarily а canonical quantiza- 
tion but, more generally, a replacement of classical quantities by operators (functions 
of position and displacement operators) then it is easily seen that the second quanti- 
zation is not the only alternative deserving the name of a field quantization. The 
other alternative is Yukawa's theory of non-local fields.? 
pU MAX AV EE rct ЭН катты т олда ma s 


! We denote by № the original Planck constant divided by 27r. 
2 Born and Peng (1943) may be regarded as predecessors. 
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This alternative type of field quantization consists in the assumption that the 
field quantities y are not functions of the classical quantities x, but (similarly as any 
physical quantity in the theory of the material point) functions of the operators x, 


and d, 
VU «= (та (6) 
Thus, the classical field quantities become quantum operators. Going over to the 
x — representation every field quantity becomes a matrix 
VE MAPAS (6) 


To distinguish this new field theory from the traditional one Yukawa introduced the 
name "non-local" field thecry. With regard to the dependence of (6' upon two 
points x' and x" it would be better perhaps to call it a "bilocal" field theory. 

(6”) is, in general, a non-diagonal matrix. A diagonal matrix would be obtained 
only if V/ were independent of the displacement operator. In this particular case we 
would have 


22 х! | y | x" Б> == р(х) д (х x. a) (7) 
and 
<x | Ч/ж | y | x" — = f < х! | yr х" = а“ x" < ж | y | х"! mm — 
= y* (x) д( — x") y(x”). (77 


This shows that the following expression may serve as a correspondence-theoretical 
transcription of densities 


ol) = fdir «x (wx (8) 
where 
рзы и (9) 


аге Yukawa's variables. 
The condition of localizability is equivalent to the following supplementary con- 
ditions : 
[0, х] = 0, и=\1...., 4. (10) 


Indeed, (10) means that У is independent of d,. 

As a matter of fact, the transition from the local to the bilocal field theory is 
a quantization. It is similar, but not identical, with the traditional canonical quanti- 
zation. Both types of quantization introduce the displacement operator but both in 
a different way. The canonical quantization introduces them by means of the condi- 
tions (5) so that physical quantities become operators 


Q = Q(x, p) where |х;, pj = thd;;. (11) 


The transition to classical theory is obtained in the limit A — 0. On the other hand, 
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the new procedure of quantization introduces into the field quantity the operator id 
on equal footing with the operator x. But these two operators have différent dimen- 
sion. Therefore we must introduce a fundamental length / and say, more ргесі- 
sely, that У (as well as any observable О being a functional of Y) depends upon both x 
апа i/?d, 

Q = ОТУ (x, il*d)] where [x,, il?d,] = И2д,,. (12) 


The transition to classical field theory is obtained in the limit I > 0. Thus, we may 
call briefly these two alternatives of quantization h— quantization, and /—quan- 
tization. In both cases the mathematical basis is the same: it is an extension of the 
classical space (x,) to the operator space (x,, id,). The traditional quantum field theory 
may be regarded as defective just because it fulfilled the program of quantization 
only partly. The two modes of quantization are not inconsistent with each other. 
We may hope that a full, and only a full, quantization of fields, being a union of the 
h — quantization and / — quantization, may provide a satisfactory theory of elemen- 
tary particles. 


3. A generalized d'Alembert operator 


So far as we are not concerned with multiplication of field quantities, we may 
regard the matrices (6’) simply as functions of two points W(x’, x") or, in terms of 
Yukawa variables (9), as functions ¥ (x, г). With the aid of (4) we may easily verify 
the following relations, valid for any F(x, а) 


<x | fd, |x" > = ЗЕ, Lu id oF |e ae 
8x, Aen 
(13). 
OF (x, г) 

д д ЕТ 

< ME Ade at = (| «x | F| "> = р 
É x 24 TELE) us 

Эти 
«xx Fx = 2, Ed c c ev 2 cx Y> (15) 

Р А, кте [err F(x, г) 
E NEST P s (16) 
We shall also use the following notations. 
| ыт Қ d) MO 

Pu = Эх, ' En SEIS NE (17) 


It cannot be a pure chance that in the local field theory every field quantity sat- 
isfies the Schródinger — Gordon equation. This equation, or the d'Alembert opera- 
tor appearing in it, must be of importance also in the / — quantized (bilocal) field - 
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theory. The question is what is the proper generalization of the d'Alembert operator 
for the bilocal field theory. In view of (13) it might be 


—4а, Фа, or 2 (а P + Уа?) (18) 


or a linear combination of both these expressions, which may be written most con- 
veniently in the form 


ЕД 


= [4 а,]_, a, + C1 [[V, 411,41, (18) 


def 

with an arbitrary dimensionless constant C. The generalized wave equation may be 
simply 

[p | - 9. (9) 


In terms of Yukawa variables (9) and with the aid of (17) this equation may be writ- 
ten also in the form 


(р> + 4C* k?) ¥ (x, г) = 0. (20) 
This equation may be solved by the method of separation of constants 
Spec а (21) 
M? 
kı = АСС; (22) 


Equation (21) is the usual Schródinger-Gordon equation in terms of the external 
variables x,. Obviously, the eigenvalues of the operators p will play the rôle of momen- 
tum and energy of the particles associated with the field. (22) is an equation for the 
internal variables r, and may be interpreted as an eigenvalue equation for the mass 
parameter M. 

In two cases: if C = 0, or if W is independent of r,, we would get the usual 
wave equation for particles with a vanishing rest mass. But a correct transition to 
the local field theory should be performed by taking account of the supplementary 
conditions (10). These conditions constitute a certain type of constraints reducing the 
number of degrees of freedom x, г,. Instead of the function ¥ (x, г) we may con- 
sider another function ф(х, г) of the independent variables only. However, this new 
function ф will not satisfy the original equation (19) or (20) but a new equation 


[p |- Re E 
where the operator R accounts for the reaction forces brought about by the con- 
straints. From (16) it is seen that (10) is equivalent to 


т, VP (x, т) = 0, ш = 1, ..., 4. (24) 
which means that (x, г) vanishes unless 7, — 0. Thus, in the limit of a local field 
the only independent variables are the х,/5, so that p = ф(х). The rôle of the term 


R is just to eliminate these of the displacement operators appearing in the generalized 
d'Alembert operator which are inconsistent with the constraints г, = 0. Such аге, 
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obviously, all the k,’s, whence R must be identical with 4 C*47 and (23) simplifies to 
the wave equation with mass zero: 


2 , 
p, 9 (x) = 0. (23^) 
The above transition to the local case is trivial, but in the next section we shall have to do: 
with another type of constraints and the procedure of eliminating subsidiary degrees 
of freedom with the help of suitable reaction forces will cease to be trivial. 


4. The supplementary conditions 


Similarly as in the case of h— quantization the time-like components of the 
operators x, and id, have to be treated on a different footing from the space-like 
components (cf. (5)), otherwise we could not find any reasonable interpretation of 
the formalism. To this end a supplementary condition 


[У, d,]., x,]_ = 0 (25) 
is needed. In terms of Yukawa variables this condition is 
T, Pu P (x, т) = 0. (25) 


For a particle at rest (р = 0) this means that У vanishes unless rọ = 0. Hence, only 
a single time variable x) = x9 = x is involved effectively in the description of the 
particle, and only space-like components of r, come into play. 

The necessity of a further supplementary condition is less obvious but may 
be made plausible by the remark that in the limit J > 0 the theory should go over 
into the traditional local one, i. e. the generalized supplementary conditions should 
go over into the conditions of localizability (10) or (24). However, by itself this require- 
ment is not sufficient to determine the form of the remaining supplementary con- 
ditions, so the question must be solved by group-theoretical considerations. Supple- 
mentary conditions mean an effective reduction of the number of dimensions of the 
operator space 0 — (x, d), since the field variables restricted by supplementary con- 
ditions depend effectively only on the operators belonging to a subspace of the eight- 
dimensional space O. Thus, the question of supplementary conditions is a matter of 
geometry: the supplementary conditions have to exhibit the intrinsic covariance 
properties of the O-space. 

As previously remarked, the "points" of this space may be subjected to Lorentz 
transformations and to canonical transformations. The supplementary conditions 
must be covariant with respect to both. The condition (25) is invariant under Lorentz 
transformations and under the group of reciprocal transformations (3) which is a nor- 
mal divisor of the group of canonical transformations*. It is to be expected that these 
invariance properties are a general feature of any supplementary condition. Their 
number should be equal to the number of (independent) Lorentz- and Born-inva- 
SS 1n E e er ee E e ыы НИ ер LOS D 


— 


*) Note added іп proofs: it is а normal divisor of the subgroup of unitary canonical trans- 
formations. ' | 


niat 
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riant expressions. It may be easily seen that there exist only two such independent 
expressions yielding condition (25) and the following one 


n [ЛА vd,,|_ qF БУ] 5 x,]- = bY, (26) 


where 0 is a dimensionless constant. By adding and subtracting (25) from (26) the 
two supplementary conditions may be written in a more elegant form (in units / = 1). 


ПР, xJ- = HS id], id] _ = b Y, (27) 
where x' and а” are given by (2). The two conditions (27) interchange under Born 
transformations. 

Since the introduction of new fundamental constants is to be avoided unless 
they prove to be absolutely necessary, we put simply b = 0. In terms of Yukawa 
variables (26) becomes (for b — 0) 

( DA 5 г) (г) = 0. (26” 

From (257) and (267 it is seen that the above correspondence requirement is satis- 

fied, and in the limit / — 0 the reciprocal supplementary conditions degenerate 
into the conditions of localizability (24) or (10). 

Condition (25’) means that, for a particle at rest, the field quantity vanishes unless 

го = 0. Then the second condition (26’) shows that it vanishes also unless r = | 7 | = 


== Үр? = l? M. Thus, the internal structure variables reduce to two angles Ф, 9 
on a sphere with radius /? M in 3-dimensional space’. The new supplementary 
conditions, being invariant (scalar) with respect to Lorentz as well as Born transfor- 
mations, are more natural than the conditions of localizability (10), and more satis- 
factory from a group-theoretical point of view. 

In the static limit the primary element in the bilocal theory restricted by the 
two supplementary conditions (25) and (26) is a directed point (the direction being 
determined by the two angles Ф, 9 on Yukawa's sphere). Therefore, the underlying 
geometry is related to that of Finsler. 


9. The rotator model of elementary particles 


Let.us investigate the solutions of the generalized Schródinger-Gordon equa- 
tion (19). The method of separation of constants yields (21) and (22). The solution of 
(21) may be Fourier analyzed*: 


1 i 2 ip их 
ЖЕ? 2 = gg) 44рд (p, + M?) Ф (p, r) еФихи. (28) 
Тһе mass eigenvalue equation (22) holds for Fourier transforms as well 
M? 
ky V(p, г) = ACA V (p, т) (29) 


ecl ИЛИНЕ cc M a e E uuu c ul RE 
3 Our second supplementary condition differs from Yukawa’s (1950) which was not reciprocally 
invariant and led to a radius r = A independent of the mass of the particle. 
4 То simplify notation we denote the operator p, and its eigenvalues by the same symbol. 
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The auxiliary conditions (25’) and (26’) may be taken explicitly into account by consid- 
ering instead of Ч (р, г) another function g(p, г) defined in the subspace of the only 
independent variables (out from the four r,). Obviously, the reduced function will 
not satisfy the original equation (29) but another equation taking account explicitly of 
the reaction forces appearing in consequence of the constraints 


M? 
|4- 46-4 ul 9 (p, г) = Ке(р,т) tus, 79, "а= 1. (30) 


The form of the term R may be easily found by remarking that these of the displace- 
ments k, must drop out which are inconsistent with the constraints. For a particle 
at rest these аге k, and k,. Writing 


о ЕЕ 
= — тн (31) 
where r = | r |, and S? is the operator of the square of the angular momentum with 
the eigenvalues 
52 —n(n--1, n - 0, 1,... (32) 
equation (30) reduces to 
S2 Mê LJ 
E е2 A 9 (p, 7) = OG =0,r=0 F1=PM] (33) 
so that 
Әз 20259459229 Se 
О EC qe sa Us 
R әй ar? pp ky — rà (34) 
for р = 0. Equation (33) yields а formula for the mass spectrum 
^» 2, 
М,- : Ұлпа ЕО п O (35) 


In the local limit M, tends to infinity except for the case n = 0. This accounts for 
the well known difficulty encountered in the local field theory. To carry out a renor- 
malization program it was necessary in the local theory to assume an infinite mechan- 
ical mass except for the case of a vanishing rest mass (e. g. photon) where both 
parts: mechanical mass and field mass were found to be indeterminate. This compar- 
ison clearly shows the singular character of the limiting case Г = 0. 

The solution of (33) are spherical harmonics S,,,(P, 0). By means of a Lorentz 


transformation we may genaralize the solution to Кы. p's. Denoting the Lorentz 
transformed spherical пай by S NI we may write 


(n) 
Фит (2, T) = Cam CPB S, „ (p) with pe ps Е — M? (36) 


In consequence of the linear character of the. wave quation and of the supple- 
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mentary conditions the general solution in the subspace of independent varia- 
bles is 


X Pan (52. (87) 
p 


The appearance of the operator S? in the eigenvalue equation (33) shows that 
we have to do here with a rotator model of the particle. The particle (in its system of 
rest) is represented by the surface of a sphere with a radius increasing linearly with 
the mass of the particle. The states of internal structure of the particle are spherical 
harmonics. The particle associated with S,,, possesses a rest mass M, given by (35). 


6. Connection with Fierz's theory of higher. spins 


Equation (33) may be written, with the aid of (34), in the form 


28 927 2 ш? 
Е Ss chop gue te roll ИНЫЕ E 


912 
Its reu for the case an arbitrary p is 


|^ ita Е E p (p, 1) = Ооа Ir (39) 


This equation is the original Schródinger-Gordon equation (20), or rather its Fou- 
rier transform, supplemented by the reaction forces brought about by the constraints 
(25) and (26). ф(р, г) is the counterpart of the function ¥ (р, г) in the subspace of the 
independent degrees of freedom. 

Let us assume that the bilocal field function ¥ is a scalar (ps scalar). Then 
also the reduced function Ф is a scalar (ps-scalar). In view of the requirement of 
Lorentz invariance the solutions of (39) must be of the form 


P(P)> т, Ф(р),» rur, 9(D)o» «°° (40) 


where the terms. Ф(р),, are tensors (ps-tensors) of rank 0, 1, 2, ... respectively. 
These tensors do not depend upon 7? since ге may be eliminated with the aid of the 
second auxiliary condition appearing in (39). The first auxiliary condition will be 


taken account of by assuming a supplementary condition 


р, 9 (p) сото = 0. (41) - 
From (40) it is seen that only symmetric tensors (ps-tensors) come into play | 
ф(р).. zo (pei t. | (42) 
Moreover, we should assume that ei the traces of these tensors vanish 
g(p) <. = 0, (43) 


(otherwise, е. g., the third term іп ye € contain a part of the same type as the 
first term). By introducing expressions (40) into equation (39) and taking account of 
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the supplementary conditions (41) — (43) it may be easily verified that they repre- 
sent solutions provided p? = — [4 r? assumes the values 


оа 262? ACE 
But this is equivalent with the previously derived formula (35) for the mass spectrum. 
The number of independent solutions is the same, so that the solutions (40) are equiv- 


alent with the solutions of equation (30). Going over from the Fourier transforms to 
the x — space we may write the general solution in the form 


ф (х, r) ms 9 (x) JE Ти qx), xl; "u r, 9 (3), sh dide (45) 
where 
1 ; 
Vires OT 1 d* pein д (p + Mn) P (puis (46) 
satisfies the (local) Schródinger-Gordon equation 
(О mp M,) 9 (x)n....u, = 0 (47) 


and the following supplementary conditions 


qx) ы ат, testo m 5 а D(a) pe = OR le) ra (48) 


The field quantities (tensors or ps-tensors) satisfying (47) and (48) are the fields 
of Fierz describing particles with integer spin values. The tensor field of rank n de- 
scribes particles with spin л. А bilocal field quantity is equivalent to an infinite set of 
local fields of Fierz with prescribed values of the masses. The /-quantization is 
а mass and spin quantization. Each of the tensor fields p(x)...,,... corresponds to 
another irreducible representation of the group of rotations. Two bilocal field quanti- 
ties: a scalar and a ps-scalar correspond to the whole group: of rotations and re- 
flections, i. e., there is a one-to-one correspondence between the members of the two 
families and the irreducible unitary representations of this group. In this way the 
bilocal fields realize the program put forward in the introduction. | 
Let us add the following remarks: a field of the form (45) restricted by (48) 


satisfies not only the complicated equation (39) but also (separately) the simpler 
equations ; 


[2 ИТ); D, k, ф(х, г) = 0 (49” 
апа ў 


[p] — 4c K?] ф(х, г) = 0 [р> = -- [4 | | (49”) 


where K? is an invariant generalization of the radial part of the Laplace operator 
1 : 
Ka к (rr kk, - 21,4). (50) 
0 


Equations (49) are not only invariant under Lorentz transformations but also under 
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dilatations. This yields an interesting connection with the conformal group. Our for- 
malism violates the conformal invariance only at this point where the ا‎ 
length J comes into play, i. e. by the second supplementary condition pie 25. 
This ш shuld be understood in (49) rather as a restriction upon d XE 
ues of p? than upon the domain of variability of ra 


7. A generalized Dirac equation 


In order to obtain a family of particles with half integer spins 1/2, 3/2, . 
may generalized the Dirac operator on similar as the d'Alembert operator and 
assume for it 


Yu Bit de s gn [5 dj. (51) 
Putting C — 1 the generalized Dirac equation takes a particularly elegant form 
y, Ud, = 0. (52) 
Its complex conjugate equation is 
d, V y, = 0, where V = V* y, (52) 
Going over to the x — representation we get, with the aid of (13), 
д , 1! 
yas x > 0), (53) 
Әх, 


ог, іп terms of Yukawa variables 
д д 
Yu 2. — 2C? z) 7 ,) = 0 or у, (ри — 262 ku) P (х, г) = 0. (54) 
2x, Әт 
where we have introduced again the constant C for generality. Going over to the 
Fourier transform V'(p, г) and taking account explicitly of the supplementary condi- 
tions (25) and (26) which are assumed to hold generally, we may investigate a new 
function (р, г) of the independent variables only which satisfies the above Dirac 
equation supplemented by reaction forces brought about by the constraints > 


Yu (b, 2 Ce к) Y (p. r) —R V (p. 1) ри =0, = 19) (55) 


The form of the operator R may be easily inferred in the rest system of reference 
(уара — 2y, C Е,) v(p.r) = RO y (p, г) 0, n, --0, п (55) 


9 ə pem | 
where the derivatives zi and 2; have to drop out from (55’) being inconsistent with 
74 r 
the constraints. 


The operator y; k; (i = 1,2,3) may be split into a radial and an angular part? 


5 Cf. Dirac quantum Mechanics (3d ed, 1947). Our J is identical with Dirac'sj while our iy, y, is 
identical with Dirac’s e.. 


rr 
are the radial component of y and the momentum cononically conjugate to r respec- 
tively. J is the well known operator 


J = в (om + 1) (58) 
commuting with y, as well as- with y, with eigenvalues + 1, + 2, ... The operator 
R must be 


Ay CE NES NES; 57 
уди ба 1 3P |- + (57) 


RO = — 2C? (y, k, + y, k,) (59) 
so that (477) reduces to 
faz + 20" i у, „г ф(р,7)--0 (60) 
or 
(һ T 2€ у, J ms) y (р, т) = 0 j-on-or-" ра (61) 


Since the eigenvalues of y, J are also + 1, + 2, ..., the eigenvalues of Гро аге + V2, 


+ үй, ... whence the mass spectrum is 


M, = S Nai k—-12,. - (62) 


On account of the appearance of the operator J, instead of S? in the mass eigen- 
value problem (60) we-may call this alternative a ,spin-rotator model“ of the par- 
ticle. 


Putting C — 1 constitutes nothing but a normalizing condition for 1. However, 
we prefer to retain this constant C which may assume different values for different 
families while / will be regarded as universal. 

With the aid of (56) the equation (60) may be written as follows 

Yapa = 2C (y; Е; — y, k,) v [р —0,r, —0,r =) (63) 
or, taking account of (57) 


„ПУ: д 
Ya Pa Y = 2C? |», ku — Y4 ky + 324 (^ rM | V ф=0т=0*= 1 (64) 


where i,j = 1,2, 3. This formula may be generalized for arbitrary р 


ТА 1 і 
Ya Pu = 263 (» ku + E Yu Pu ky Py — 3 Yn Ta Ts hy ы Уи ra) V iruPu=0 r, 7 — pd. 
б | 


% e 
| (65) 
The equation (65) is nothing else but the Dirac equation (52) supplemented by re- 
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action forces arising from the supplementary conditions (25^) and (26) while v is 
the counterpart of V/ in the subspace of independent variables. The equation (65) 
may be solved by the method of separation of constants 


(РЕМИ (66”) 
7 О кроты) — iM * 
7 u - и М» Py ply Fy 7 Уз а” (66) 


(66°) is the usual Dirac equation with the mass parameter M while (66) is the eigen- 
equation for the mass parameter. The eigensolutions of (66") are again 


vip), Ти vip), D, r, (р), РЫ (67) 


provided the spintensors appearing іп (67) are symmetrical in every pair of tensor 
indices and satisfy the auxiliary conditions 


Pa (PJ  — 0, y, wp), ++ = 0. (68) 
By introducing the terms (67) into (66”) we find again the mass eigenvalues (62) 
where the indices Ё = 1,2, ... refer to the separate terms (67) successively. The 


particles with mass M possess spin / = k — 1/2. Denoting the mass terms by the 
spin index j we have 


C s~. 1 3 


Going over to the x—space we may write the general solution of the Dirac equation 
(52) or (54) supplemented (or rather reduced) by the auxiliary conditions (25) and 
(26) in the form 


p(x, r) 3 v (x) T T, v(x), TE 7, Т, y (x) w ЕЗ MU (70) 


with the supplementary conditions 


ADAC) rige SD, 52% (е)... = 0, у, (2) = 0. (71) 


The bilocal field quantity restricted by the reciprocal supplementary conditions is 
equivalent to an infinite series of local fields with half integer spin values. 


8. Alternative types of field equations 


It is well known that particles with integer spin may be described not only by the 
Schródinger-Gordon equation but alternatively by the first order Kemmer equation, 
these two modes of description being equivalent within the local field theory. The 
question is whether the two alternatives remain still equivalent within the bilocal 
field theory We shall show that the answer is negative. The family of particles described 
by a generalized Kemmer equation leads to a mass spectrum of the same type as (62) 
derived from a generalized Dirac equation. On the other hand, as pointed out by 
Cap, (1954) it is also possible to describe particles with half integer spins by means of 
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second order equations and this description is equivalent to the usual one within 
the framework of the local formalism. We shall show that also this equivalence 
breaks down in the bilocal field theory. Let us discuss briefly these alternative cases. 

The Kemmer equation differs from the Dirac equation merely by the replacement 
of gamma matrices by Duffin-Kemmer beta matrices. In close analogy with (66) we 


may assume immediately 


(iB, p, + M)g = 0 (72^) 


1 А 7 
B |^ + го (14 р, k, p, m" k, mm Vi ф = з Ф [p ru = 0р'и = — liru]. (7251 
Го 


Тһеѕе equations may be solved in the form 


9x. r) = фа(х) + T, Ф.х), T T, т, Pal*) wv Hs ps (73) 


The proof that (73) is a general solution of (72), provided the separate solutions are 
subjected to the supplementary conditions 


Ф (2x)... vem UX) unge O ф(х). = 0, В, (о... = 0 (74) 


and correspond to the mass eigenvalues (62) or (69) respectively, is exactly the same 
as in the case of the Dirac equation and does not involve the commutation properties 
but only the eigenvalues of the gamma or beta matrices. Excluding the eigenvalue 
zero of the beta matrices®, we are left with the mass spectrum of the type (62). 
Using five row and column representation of the beta matrices ф(х) possesses five 
components necessary for the description of particles with spin zero. The next 9. (x), 
restricted by the conditions (74) possesses ten components which is necessary for the 
description of particles with spin 1, and so on. Thus the spin of the particle with mass 
M, given by (62) is n=k-1, and we may write 


M, — © ү) n = 0,1, ... (75) 


where n. denotes the spin of the particle in question. Formula (75) differs essentially 
from (35), whence the two modes of description of bosons are not equivalent. 

А possibility of a theory of particles with half integer spins described by a second 
order wave equation may be best seen by taking advantage of the Van der Waerden 
spinor notation (Bade and Jehle 1953). Let us assume (Cap 1954) 


(DI my = 9 (EIE үрк; (76) 


where y, is a two component spinor u = 1,2, и = 1,2. Taking (76) as a starting point 
for the bilocal generalization of the theory, we shall obviously obtain a family of 
ferniions with a mass spectrum of the type (35) with the only difference that the mass 


* The field quantities belonging to the value zero are dynamically dependent. 
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M, corresponds now to a particle with spin j = n + 5 . Labelling the mass eigenvalues 


by the spin index we get in this case 


1255 
М; = = V(2j)? – 1 j= = ә әсе (77) 
It may be noticed that the same result may be obtained also in а different way. 
We may start with the generalized Schródinger-Gordon equation (39) which may 
be split into two parts 


(p; + M?) e (p,r) = 0 (18) 

М \2 
Е а= E (14 (ри ky)? са [Л k, a 21 Wy ky) — [ва | ф (p, r) = 0 [Puru —0.7*,, — М3, 
(78") 


(78’) is the usual wave equation while (78") is the mass eigenequation. Now, we may 
linearize (78) 

(i y, p, + M) g(p.r) = 0 (79) 
while leaving, ex definitione, the mass eigen-equation as it stands. This partly linear- 
ized wave equation obviously leads to a family of particles with half integer spins but 
with a mass spectrum of the same type (35) as that previously obtained for bosons 
satisfying a generalized Schródinger-Gordon equation. 

Concluding it may be stated that there exist (at least) two types of mass spectra 
for bosons: (35) and (75), as well as two types of mass spectra for fermions (69) and (77). 
The mass spectra involving square roots belong to families described by the first 
order wave equations, while those with fourth roots belong to families described 
essentially by second order wave equations. 


9. The meaning of the h- and l-quantizations 


We have shown that a bilocal field is equivalent to an infinite set of local fields 
with higher spins and masses. This result is most gratifying. It shows that all irreducible 
unitary representations of the group of rotations have been taken into account on 
equal footing, in contradistinction to the traditional field theory which considered 
a limited number of fields corresponding to some arbitrarily chosen representations. 

The equivalence of the bilocal (J-quantized) field to a set of local field solves, 
in principle, the question of observables, conservation laws, and second quantization 
(h-quantization) by analogy with the corresponding problems within the local field 
theory. The separate local fields are orthogonal to each other’, so that observables 
like energy, momentum, etc... for free fields split into an infinite series of correspon- 
ding quantities for the separate local fields. These fields may be subjected successively 


1 This follows from a théorem about functions corresponding to different irreducible representations, 
or different rows of irreducible representations. 
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to the process of second quantization with a subsequent limit transition to the case 


of infinitely many fields. 

In view of the fact that spins and masses increase to infinity within each family, 
the formalism offers a method of approximations for solving problems of interaction. 
In states with not too large total energy and angular momentum only a few particles 
play an important róle. Іп practical computations we may limit ourselves to the first 
few particle sorts. Taking into account even a few of them, a considerable progress 
may be achieved in the problem of convergence: there exist several indications that 
we shall have to do with mutual compensations of infinite self energy terms (realistic 
regularization). 

The results of the foregoing sections show clearly the rôle played by the /-quanti- 
zation as compared with the h-quantization. The (first) h-quantization is, in principle, 
a quantization of the dynamical states of a single (individual) particle. Of course, it 
may be extended to more than one particle (by means of an extension of the number 
of coordinates, forming a configuration space) but this extension does not work 
automatically, it has to be performed. Then the number of particles is always definite 
and cannot change. 

On the other hand, the second quantization (quantum theory of local fields) 
is, in principle, a quantization of the number of particles (of a given sort). Now, a field 
quantity describes (instead of a single particle) any set of particles in arbitrary dynam- 
ical states. The number of particles may have or not a definite value (the state may 
be or may not be an eigenstate of the number of particles.). Thus, the -quantized 
field theory is no more а theory of individual particles but a collective theory of a whole 
sort of particles. Of course, it may be extended as well to more than one particle sort 
(by introducing new components of the field, i.e. new field quantities) but again this 
transition is not automatic, and the number of field quantities and particle sorts is 
fixed. 


Given any number of particles, the states corresponding to a different number 
of particles are (in the framework of the first quantization) simply ignored, while in 
the quantum field theory they cannot be ignored but contribute (at least virtually) 
to the developement of a given state in the course of time. In the framework of the 
second h-quantization not particles but particle sorts are treated individually; given 
any set of particle sorts, the states involving other particle sorts are simply ignored. 


Then comes the /-quantization being a quantization of the sorts of particles. 
In this case even the particle sorts are no more treated individually but the whole 
family is described by a single bilocal field quantity. States involving any other particle 
sorts belonging to the family cannot be ignored but contribute (at least virtually) to 
the development of a given state in the course of time. 

We have to do here with successive stages of a progress towards more and more 
general formulations at the cost of avoiding individual (and accidental) distinctions. 
In view of the above arguments we hope that, although some details will have to 


Introduction to the Theory of Elementary Particles 555 


undergo a revision, the general direction of progress sketched in this paper is right 
and should be useful for a qualitative and quantitative explanation of several problems 
connected with elementary particles. 


CHAPTER II. COMPARISON WITH EXPERIMENT 
l. General remarks 


Let us check the ideas concerning elementary particles developed in Chapter I 
by the present experimental evidence. First of all we shall try to systematize quali- 
tatively the known particles into families, secondly we shall compare the formulae 
for the mass spectra with the present experimental data. Of course, we cannot expect 
an exact numerical agreement between the theoretical and the experimental values of 
the masses since the former are merely mechanical masses, while the latter include 
also field masses originated from the interaction with other fields. So long as a satis- 
factory theory of interaction is lacking, the self masses (causing so to say a fine structure 
of the mass spectra) may be introduced only in the form of phenomenological correc- 
tions. The mechanical masses appearing in the /-quantized field theory depend upon 
the ratio С/1. The ratios of the masses depend merely on the ratios of the respective 
constants C. In principle we may dispose freely of this constant, i.e. assume as many 
different C's as there are different existing families. However, these constants should 
be connected in some way or other, and the formalism will be convincing only if we 
shall have a small number of different C's related by a simple algorithm to each other. 
In particular, the formalism of Fierz makes it plausible that bosons and fermions may 
be collected into ,,super-families“ with spins :0, !/5, 1, 3/ , ... so that we may expect 
the same values of C for fermions as well as for bosons. In fact, we shall show that only 
two values of C with a simple relation 1 : 2 will do to account for the main features of 
the mass spectra of the known particles. 

Now, the problem is as follows: how many particle families are expected to 
exist in Nature and what are their distinctions. First of all we may distinguish between 
families of particles according to parity and then according to whether they are described 
by a first or a second order wave equation. Thus, we may expect the existence of four 


families of bosons: 


tensor ps-tensor 


уи 6 
XN | 


first order eq. second order eq. 


apart from some exceptional cases (photon, graviton) which seem to belong to the 
value C — 0 and have to be treated separately. 
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In the case of fermions, we have again the two alternatives of first and second 
order wave equations but a more complicated situation concerning parities. As pointed 
out by Racah (1937), and discussed in detail by Yang and Tiomno (1952), there exist 
two objectively different kinds of spinors с under inversion of space axes 
according to 


y = ya (а) yp = iyw. (В) 


For a single field we cannot distinguish between a plus or a minus sign since spinors 
are two-valued functions, but we may speak about relative parities of two fields whose 
product yields either a plus sign under an inversion in case (a) (and a minus sign in 
case (В) or vice-versa. We shall denote these alternatives by (+), (+), (+i), and 
(Fi). The charge conjugate of a field quantity of type (+) has an inverse relative 
parity (F) while the charge conjugate of a field quantity of type (+t) has the same 
parity (Fi). Therefore, we should expect to have four families of fermions of type 
(+i), according to the scheme 


кк (Һ) 
ЙЫ Y 


first order eq. second order eq. 
but only two families of type (+), according to the scheme 


(+) 
x^ S 


first order eq. second order eq. (c) 


since the opposite two cases with (F) are already realized by the respective antipar- 
ticles. Curiously enough, we shall need for the cases (b) and (c) two constants c 
differing by a factor two, namely C, : C, — 1:2. This must be somehow connected 
with the fact that the group of transformations in case (а) is isomorphic with (1), 
(—1), while the group of transformations in case (В) is only homomorphic with (1), 
(—1), and isomorphic with the group (1), (7), (—1), (—i). The inverse C7" is just pro- 
portional to the number of elements of the group governing the transformation law 
of the field quantity under inversion, 


2. Boson families 


The recent analysis of the decay data of the heavy meson tau shows that most 
probably it is a pseudoscalar (Vitale, private communication). Thus, we have two 
different ps-scalar mesons pi and tau, and therefore we may expect them to belong to 
two different types of field equations. As tau is much heavier than pi, we suppose 
that pi is the lowest member of a family satisfying the Schródinger-Gordon wave 
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equation leading to the mass spectrum (35), while tau is the lowest member of the 
other ps-scalar family satisfying the Kemmer equation and leading to the mass spectrum 
(75). If it is really so, the mechanical mass of the pion is zero and its whole mass is to 
be regarded as a field mass. This seems plausible for the following reason: the ratio 
of the electromagnetic mass to the self mass due to nuclear forces (273--264)/264. 
agrees roughly with the ratio of the respective coupling constants e?/g?. We cannot 
compute the masses of higher members of this family without knowledge of the con- 
stant C/l but, at the moment, it is of little interest since the higher members of the 
pion family are expected to decay into three pions іп an extremly short time ~ 10-22 
sec and escape direct observation owing to a large value of the nuclear forces coupling 
constant. 

Going over to the tau-meson family we may abandon, in the first approximation, 
the self mass and regard its total mass as identical with the mechanical mass. Under 
this assumption we may compute C/I from (75) by putting M, = 965. The next masses 
become M, = V 2M, = 1360, М, = 1670, etc... Unfortunately no sufficiently reliable 
experimental data are available to prove or disprove these predictions. 

The mesons 0? апа X^ are most probably scalar. The striking similarity of the 
value of their masses with that of the t-meson is easily explicable within our scheme 
by assuming that they constitute the lowest members of a family described by scalar 
and ps-scalar bilocal field quantities and satisfying the Kemmer equation. with the 
same value of the constant C. Their long life-times may be explained by. assuming 
that they are coupled to the nucleon field by a coupling constant g by several orders 
of magnitude smaller than the coupling constant of nuclear forces (coupling between 
nucleons and the pion family). To the question of their copious production we shall 
come later. 

The results are put together in the following table. 


Table 1 
ps-scalar scalar 
I-st order eq. I-st order eq. 
spin 0X1 spin QN IG S 
parity RR party J 
mass (965) 1360 1670 mass 965 1360 1670 
name әйтте омыр? name pen ? 
ps-scalar scalar 
IL-order eq. II-nd order eq. 
spin 0 1 2 spin 0 1 2 
parity rah pane parity + — + 
0 
md (270) ? е mass 0? 


name xr (short lived) name Pal, ? 
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The values of the masses written in paranthesis are adjusted, the remaining ones 
follow from the formalism. We notice the lack of any evidence of particles belonging 
to the second order equation and scalar. 

The arguments in favour of our systematization of bosons аге: (7) Experimental 
evidence yields two ps-scalar mesons: pi and tau, one of which is much heavier than 
the other in agreement with the theoretical prediction. (ii) Experimental evidence 
yields scalar particles theta with about the same mass as tau which fits well into our 
scheme without necessity of introducing any new value for the constant C into our 
formalism. (iii) The observations of Leighton? of the particle V3 = 7) with the probable 
decay scheme 7° ->л + K + 60 Mev. 


3. Light fermion families 


Our first task in systematizing the fermions is to look for a particle possessing 
a mass close to that of tau. If such a particle exists, it should possess spin !/; and belong 
together with т and @ to a super-family of particles with spins 0, 1/;, 1,9/5,... described 
by a first order equation with a common constant C. We suspect that it is identical 
with the particle K with a mass 935 + 20 and most probably with a two-body decay. 
If our supposition of К+ being a fermion is right, the neutral decay product is a light 
boson К} —* + 1 n9. The idea that n? is a ,scalaron*, i.e. the lowest member of the 
missing PUR family of bosons seems to us very attractive. The difference between 
the mass of the tau meson 965 and the meson K, may be explained by an additional 
self-mass in the case of bosons. The meson К+ exhibits a similarity with the proton 
in that it is always positively charged. Another striking fact is that its mass is nearly 
one half of the protonic mass. This may be due to the circumstance that the constant 
C is twice as large for heavy fermion families than for light fermion families (cf. Sec- 
tion 1). We may-arbitrarily put C = 2 for the nucleon family which merely normalizes 
the constant J. Applying our formula (69) with C = 2 to the proton (M = 1936), 
we get a reasonable value for the universal length J = 0,594. 10-13 cm. Then the light 
fermions are characterized by the constant C — 1 and the theoretical value of the 
mass of the meson K7 becomes 918. The higher mass values for this family follow from 
(69): 1300, 1590, etc... 

Besides the К+ family characterized by C = 1 and a first order equation, we 
have the other alternative of a second order wave equation with the mass spectrum (77). 
This yields M — 0, M — 1090, M — 1430, ... Obviously the value zero will apply 
to the lightest fermions: electron, neutrino, muon, while their relatives with higher 
spins are expected to have masses beyond one thousand electron masses. 

Let us go over to the question of systematization with regard to reflection symmetry 
The analysis of Fermi-interactions given by Yang and Tiomno led to the 
conclusion that most probably the nucleon is characterized by a (+) parity while 
light fermions are characterized by a (+i) parity. Assuming this viewpoint, we have 


8 Cf. the article of Hodgson (1953). 
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to expect the existence of four families of light fermions altogether: two families of 
fermions described by first order wave equations with reversed relative parities (+i) 
and (+i), and two families of fermions described by second order wave equations with 
reversed relative parities. In particular we sheuld expect, besides К}, another particle 
with spin !/; and about the same mass, and two types of fermions with spins !/, and 
masses. close to the value zero. As for the lightest fermions we have three instead of 
two kinds of particles. This may be understood in.the following way: the neutrino is 
simply a neutral counterpart either of the electron ог the и-тезоп. The small mass 
difference between the electron and the neutrino favours the first alternative, however, 
since the mechanical mass of the j-meson is also zero, the above argument is not 
very strong. Àn analysis of the interactions of the lightest fermions with other particles 
favours rather the opposite alternative: the neutrino seems to have the same parity as 
the u-meson and opposite to that of the electron so that у and и belong to the same 
family labelled e.g. by (+i) while the electron belongs to the other family labelled by 
(Fi). To see this let us consider the pion decay 
nE ج‎ р + >. 

If the parities of u and у are both (Fi) then the product of both parities is minus 
(Fi) (Fi) = (—), ie. the same as the parity of the pion. If we assume a reversed 
parity (+7) for the electron, then the alternative decay scheme л= — e+ + v (which 
has never been observed) is characterized by a difference of signs on both sides. Of 
course, this is not an absolute selection rule since the conservation of parity may be 
secured by assuming that the decay products carry off an angular momentum h and 
form a P-state?. Nevertheless, the decay without carrying off an orbital angular momen- 
tum seems more probable and so the decay into и + v is favoured in comparison with 
the decay into e + v. 

The same rule holds for the ш-сарішге favouring Р? + u^ — N + v in comparison 
with Pt + u^ Р" + ег, which would be possible only if В+ + ег formed а P-state. 

We may expect that the values of the masses of the charged members of the 
u-family will be similarly shifted (roughly by about 200) as is the case with the muon 
itself, so that we may expect that the charged counterpart of the muon with spin 3/2 
possesses a mass M~ 1300 instead of 1090, the next MS? — 1600 instead of 1430 
and so on. The same rule of conservation of eigen-parities (the decay in an S-state 
without carrying off an orbital angular momentum) favours the decay schemes 


К+ > и= фу + > 
for the particle with spin 3/, belonging to the muon family, and 
К+ >et yp} ty = 
for the decay of the corresponding particle belonging to the electron family. Such 


9 The law of conservation of angular momentum is not invalidated since both spins may form а ге- 
sulting angular momentum opposite to the orbital angular momentum. 
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particles with masses beyond 1000 have been actually observed and denoted by K 7 
(ог x*) and К+. Their masses have been estimated аз 1000—1300 (Glasgow Conf., 
July 1954)19, 


The assignments of light fermions are given in Table 2. 


Table 2 


Ci 
I-st order eq. 


(+i)-family (Fi)-family 
spin 1 3/5 Р spin 1/5 3/5 5/5 
parity аі - i parity Fi +i F 
mass (918) 1300 1590 mass 918 1300 1590 
name ҚТА? D ? te 4 


II-nd order eq. 


spin 1/5 3/5 2/5 spin 1} MP 5, 

parity +i E +i parity Fi +1 SE 
0 1090 1430 

or, tent (207) 089857591300 8591600 mass 0 1090 1430 
y ж 

name ut Ae ? name е Ke ? 


No particle belonging to the (-Fi)-family described by the first order equation has been 
observed. It may be expected that the lowest member of this family should decay 
into an electron and a neutral boson. The results of our analysis enables us to classify 


not only the lightest fermions but also the remaining particles discovered recently 
in the mass interval about 1000 m,. 


4. Heavy fermion families 


As pointed out in Section 1, we have to expect only two further families of fermions 
of the type (+) governed either by the first or the second order wave equation since 
the alternative (F) is realized just by the respective antiparticles. The particles belong- 
ing to the (+)-families cannot decay into those belonging to the (2-)-families. This 
is an absolute selection rule: the law of conservation of the number of particles minus 
antiparticles for the two types of families separately. 

We assume that the nucleon is the lowest member of a family described by the 
first order wave equation with the constant C = 2. The masses of the higher members 


10 Photometric measurement of K+ yields 1050+100 in good agreement with our 1090 (Friedlander, 
1954). 
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of this family are given by formula (69) and, of course, are twice as large as the masses 
of the K} -family. М, = 1836, М, = ¥2. 1836 = 2597, Му, = V3 1836 = 3180, 
etc. The value 2597 agrees very well with the experimental value 2582 + 10 for the 
cascade hyperon Y. Moreover, we notice the following interesting fact: since the 
parities alternate within each family, the hyperon with spin 3/, and mass about 2600 
possesses a reversed parity to that of the proton, i.e. the same parity as the antiproton. 
On the other hand, the sign of the charge of the cascade hyperon is the same as 
that of the antiproton. It is hard to believe that this coincidence is merely accidental. 

The masses of the other family of fermions governed by the second order wave 
equation are given by formula (77). Assuming the same value C — 2 as for the nucleon 
family we get in this case 


4 
М, = 0, My, = VZ. 1836 = 2183, М, = 2860 еіс... 


The value M., agrees surprisingly well with the experimental value 2181 for the 
hyperon 4°. It is hard to believe that the coincidence of the theoretical values 2597 
and 2183 with the experimental values 2582 + 10 and 2181 + 1 is only accidental. 
These values have been obtained without any adjustment of constants except for the 
simplest assumption that the ratio C/I is the same in both cases. We are inclined to 
regard the above numerical agreement as a strong confirmation of the principal 
correctnes of the ideas underlying our formulae for the mass spectra. 

Of course, the formalism in its present form is not yet completely satisfactory. 
In order to prevent the decay of nucleons and hyperons into the light fermion with 
mass zero appearing in the A°-family, we have to exclude the existence of this par- 
ticle by an additional rule: The index j should run in this case from 3/, and not from !/,. 
This is an additional rule beyond the formalism itself. Alternatively, it is conceivable 
that the future theory of the coupling constants may prove that the coupling „соп- 
stant“ is a function of the spin, indices, assuming a value zero for the value j = 1); 
of the index of the A°-family, preventing any interaction with this particle. 

The light particle being excluded in this way or other from the A°-family, we may 
go over into details of the interaction involving nucleons and hyperons. Assuming 
that the parity of the hyperon 4° is the same as that of the nucleon the decay scheme 


A? Р + л 
is consistent with the combined laws of conservation of angular momentum and parity. 
The decay products must unavoidabaly form a P-state since they have to carry off an 
orbital angular momentum A. Since Y- has a reversed parity to that of the nucleon 
the decay scheme 

Ү > Л° Hw 
is also consistent with the laws of conservation of angular momentum and parity if 
the decay products form an S-state. On the other hand, the decay scheme Y > N+ m- 
(missing in experiment) necessitates the decay products to form a D-state to preserve 
the conservation laws. But the necessity of carrying off an orbital angular momentum 
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2h means rather a strong selection rule, which seems to suffice for the explanation 
of the cascade character of this hyperon. Thus, the assignment of the cascade hyperon Y^ 
as the particle with spin 3/, belonging to the nucleon family and 4° as belonging to 
the other family finds a confirmation not only in the close agreement of their masses 
with those predicted theoretically but also in the analysis of the possible decay schemes. 

We notice that the value of about 2320 for the charged hyperon Q* is not found 
among the theoretical mass values. The decay scheme of this hyperon is similar to 
that of 4°, so that it seems natural to regard it as a charged counterpart of A°in 
spite of an unexpectedly large difference between their masses. This difference is 
analogous to that between the muon and the neutrino. We are unable to explain why 
sometimes the electromagnetic masses should be so large. It is not excluded that they 
come about from higher order terms involving the interaction with pions. Such terms 
could accidentally, so to say, compensate in the case of proton by a fortunate difference 
of signs in some self-energy terms. If it is so, then we should expect that also the mass 
of the charged higher member of the A°-family will be shifted from 2860 to about 
3000. These two values seem to find confirmation by experiment. We would like to 
call attention to the event of Fry and Swami (1954) who observed a decay of a nuclear 
fragment with the emission of a pion with an unusual release of energy Q — 360 MeV. 
This corresponds to the mass value of about 2820 for the hyperon (let us call it F) 
in agreemerit with our value 2860 for the heavier hyperon belonging to the A°-family. 
On the other hand, Friedländer (1954b) observed an event which may be interpreted 
most naturally as a decay of a charged hyperon into a pion and a neuiral particle. If 
this latter is a neutron, then the mass of the primary is about 2950. We notice that 
the theoretically predicted masses show a systematic shift towards larger values as 
compared with experiment. The results concerning fermions of the (+) type are put 
together in the following table. 


Table 3 
C2 


Nucleon family 1-51 order eq. A°-family П-па order eq. 


Чұ 


х 2/5 5/3 
parity dE Er + 
(1838) 2185 2860 
mass (1836) (2320) ~ 3000 
0 .0 
N,P + 4 
name Qt FX 


9. Some final remarks: 


We should like to stress the following interesting fact: the law of conservation 
of the number of particles minus antiparticles shows that Q- cannot be an antipar- 
ticle of Q*; (being fermions they must have antiparticles of their own). This shows 
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that Q+, Q- and 4° form a triplet corresponding to the value 1 of the isotopic spin. 
Thus, instead of the old rule that the first (second) order field equations describe par- 
ticles with half-integer (integer) spin (which has been disproved) it seems that there 
exists an analogous rule for isotopic spins. The nucleon possesses isotopic spin 1/, 
while the hyperon A possesses isotopic spin 1. Now, the question arises whether 
the same rule holds for fermions of type (i)? The meson K, is described by the first 
order wave equation and is always positively charged. Assuming the existence of a neu- 
tral counterpart its isotopic spin could be !/,. The electron stands alone, so that its 
isotopic spin is zero in agreement with the pretended rule. On the other hand, the 
rule would hold for u and v only if there existed two different muons or neutrinos. 
If this rule holds true (generally) it may contribute to a future understanding 
of the notion ,,isotopic spin“, whose meaning is still so obscure nowadays. 
Finall let us say a few words about the question of the long life times 
and copious production of heavy particles. As pointed out by Pais (1952), the 
apparent contradiction between the long life times of heavy bosons and hyperons and 
their copious production may be resolved by a *two-coupling-constants hypothesis*. 
We have to assume (at least) two coupling constants: a large G and a small g in the 


case of couplings of the type шурф, where p, denotes а boson field quantity while 
v; denotes a field quantity describing nucleons or hyperons. We notice the following 
interesting reciprocity: if o, denotes the pion then the coupling constant is С for 
| = j (or g for i Æj, that is for different spinor particles). If, on the other hand, 
9, denotes O° (and probably also т) then the coupling constant is g for i Aj (or С 
for i j). This rule accounts not only for the long life times of heavy bosons and 
hyperons but also for a simultaneous production of pairs of heavy particles (a heavy 
meson produced simultaneously with a hyperon, or two hyperons produced simulta- 
neously). Moreover, this rule accounts also for the share of hyperons in nuclear forces: 
he A-nuclei discovered by Danysz et al. 


KPATKOE СОДЕРЖАНИЕ 


Райский, Введение в теорию элементарных частиц 


à 


`Обобщенный полевой формализм (тесно связанный с формализмом Юкавы) 
интерпретирующий фундаментальную длину 1 можно рассматривать как нового 
рода квантизацию — „квантизацию 1” — в отличии от классической „квантиза- 
ции |”. Оказывается, что квантизация | является квантизацией массы и спина 
и ведет к вращательной (или спино-вращательной) модели элементарных частиц. 
Частица представлена сферою (собственно поверхностью ее) с радиусом ЁМ (где 
М — масса покоя). Общее решение управнения поля (первого или второго по- 
рядка) включает целое семейство с внешними спинами и массами. В работе рас- 
сматривается близкая связь с теорией Фирца высших спинов. Частицы класси- 
фицированны в семействах и сверх-семействах соответственно типу уравнения 
поля (1-го или П-го порядка) и типу четности. Из общих соображений согласно 
с теорией групп следуют четыре семейства бозонов и шесть семейств фермионов, 
Существует подобие формы между всеми однозначными или многозначными не- 
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приводимыми представлениями групп вращений и инверсий и возможными типа- 
ми элементрных частиц. 

Все роды частиц известных в настояшее время (легкие фермионы, легкие 
и тяжелые мезоны, нуклеоны, гипероны) подходят к нашей схеме, с исключе- 
нием фотонов (и гравитонов), которые необходимо рассмотреть отдельно. Каче- 
ственные свойства (схемы распада и т. д.) можно выражать посредством правил 
отбора исходящих из сохранения момента количества движения и четности. 
Теоретические значения масс, которые следуют из этого формализма, в большин- 
стве случаев замечательно хорошо согласны с экспериментальными данными. 
В остальных немногих случаях (м, л, QE) согласие только грубое, но это pac- 
хождение возможно вызвано „тонкой структурой” уровней массы (собственных 
масс) обусловленной взаимодействием. 
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ON THE SURFACE TENSION OF IDEAL SOLUTIONS 
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(Received January 24, 1955) 


In a recent paper Puchalik (1954) has reported an experimental study on the 
surface tension of some solutions, made in order to investigate the association processes 
accompanying dissolving of one substance in another. It seems, however, that the 
surface tension is influenced by association processes in a very complicated manner 
and before drawing any conclusions from such a study changes in the surface tension 
of ideal solutions must be first determined. 

For ideal solutions Freundlich (1922) derived an equation based on the additivity 
of surface tension of pure components. Stachorski’s (1928) empirical equation for 
a solution of two components is 


0,03 
a(l — а) + оз’ 0) 
` where 7, сі, сз are respectively the surface tensions of the solution and of its first and 
second components and c, is the molar concentration of the first component. Denoting 
the molar concentration of the second component by c; = 1 — ды we can easily trans- 
form (1) into 


С = 


Бае (2) 


which shows that equation (1) is based on the penc of reciprocals of the surface 


tension. 
Hammic and Andrew (1929) have assumed that the parachor of an ideal solution 
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can be calculated by taking into account the additivity of this constant: 
Р = aP, + CP, (3) 


They have determined the parachor P of the solution by measuring the surface tension 
and the density of a solution of known concentration and comparing P with the value 
of the parachor calculated from equation (3). The same method was adopted by Pucha- 
lik. Thus for determining the parachor of the solution two physical constants have 
been used, both of which can be influenced by association: the surface tension and the 
density. It is well known that in non-ideal solutions the density depends considerably 
on concentration. Because of this fact the method of Hammic and Andrew does not 
allow us to determine the influence of association on the surface tension. To investigate 
this influence we have to calculate the surface tension of the solution from the physical 
constants of pure components only and compare it with the measured value. 

An appropriate equation has been obtained some years ago by the author in an 
d work (Mierzecki 1949)* starting from equation (3). Taking P = Мо Ча 
= И 010 ^(M — molecular weight), we have for an ideal solution of two components 


1 
оза ат e (1 б^, то1, AR 957 Же (4) 
For an ideal solution there is по contraction and therefore . 
à M М, 
Ж к= Vols i Co Жа == Сү » «ir Сә = 
Finally 
M, М.А 
co,” Em d C024 2 
g = | — غ‎ (5) 
M, 


M. 
e 


represents the normal behaviour of o as function of the concentration c. 
Puchalik assumes such a function to be given by Szyszkowski's equation (1908) : 


o = og — b In (1 + ас) which is valid for aqueous solutions of superficially active 


substances and contains a and b, two constants different for different substances 
dissolved. Such solutions can in no way be considered as ideal solutions and 
cannot therefore serve as basis to draw conclusions regarding association processes. 

From equation (5) we see also that one cannot assume that for a given solution 
the dependence of c on c is linear and at the same time the parachor is an additive 
function of the concentrations. Nevertheless Puchalik writes: ,,Die Kurve für die 
erstgenannte Lösung (i. e. for the С,Н,ОН — C,H, solution) lässt sich gut durch eine 
gerade Linie approximieren und bei der gegenwürtigen Genauigkeit der Messungen 
kann man nicht bestimmt entscheiden, ob die vorhandenen Abweichungen von der 
. Linearitát wirklich bestehen“ and at the end of his paper we read: „Nur im Falle der 


А 
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Aethylalkohol — Benzol Lósungen erweisen sich die Parachore als streng additive 
Funktionen der Komponenten.* 

It seems that the exactitude of the contemporary measurement technique of sur- 
face tension does not allow Puchalik to determine in an definitive manner either the 
linearity of the о = (с) function or the additivity of the parachor. To day it seems 
also impossible to find out-experimentally which of the equations is true: that of 
Freundlich, equation (2) or equation (5). The values of surface tension of the compo- 
nents of ideal solutions cannot differ strongly and in that case the experimental errors 
are greater than the differences between the values calculated from these equations. 

The problem of the influence of association on surface tension is yet far from 
beeing solved and it seems that Puchalik is right in stating: „Die endgültige Entschei- 
dung ob man die gewonnenen Resultate zur Unterscheidung der verschiedenen Asso- 
ziationsarten verwerten kann bleibt künftigen Untersuchungen vorbehalten." 

The author wishes to express his gratitude to Professor A. Dorabialska for 
permission of publishing some results of his dissertation. 
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ON THE MEASUREMENTS OF 
COSMIC RADIATION IN HORIZONTAL DIRECTION 


Ву М. Kocur AND А, OLEŚ 
General Physics Laboratory of the Mining Academy, Kraków 
(Received March 4, 1955) 


There are in the literature very few data concerning the measurements of cosmic 
radiation in horizontal directions. Clay (1949) found in horizontal measurements 
a great excess of photons over ionizing particles; from the data given by Clay one might 
conclude that they refer to photons of low energy (< 10 MeV). Rogoziński (1948) 
stated that the horizontal cosmic radiation should be very penetrating. He found that 
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Scale 1:10 


Fig. la 


this radiation is practically not absorbed even by a layer of 15 cm of lead, but going 
through an absorber it produces a secondary radiation which is strongly absorbed, 
like high energy electrons. However, in both above experiments the telescopes were not 
sufficiently protected by anticoincidence counters against the influence of radiation 
coming down from other directions. The aim of this letter is to mention some results 
of an experiment which has been carried out by means of a horizontal telescope 
thoroughly protected by anticoincidence counters against different side effects. 
The counter-apparatus used is shown in Fig. la, b. The telescope registered the 
coincidences from counters 1, 2, 3. It was protected with the anticoincidence counters 4. 
To evaluate the non-ionizig horizontal component we have measured the transition 


(369) 
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curves in the lead absorber S of variable thickness ranging from 0.025 cm to 5 cm in 
the presence of the additional anticoincidence counters 5. Another lead absorber P; 
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Fig. 2 


10 cm thick, was used to eliminate the background due to photons which do not fall 
from horizontal directions. 
| The transition curve obtained with a telescope of brass-counters for the non- 
-onizing component is given in Fig. 2 (curve c); it resulted by subtracting curve b 
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from curve a. Curve a gives the dependence of the rate of anticoincidences (1, 2, 3) — 
(4, 5) on the thickness of absorber S for measurements without’ the absorber P. Curve 
b gives the same dependence for measurements with the lead absorber P. The above 
anticoincidences may be due to the inefficiency of counters 4 for registering the 
particles produced in R and absorbed in S. These particles might be decay-electrons 
of the u-mesons which arrive at В without operating the anti-counters 4 or knock-on 
electrons of the same origin. 


Ал | 
The transition curves obtained by: 


Jánossy, Rossi 
کے کے‎ Trumpy 
=o authors with brass counters 
—Ó 6v authors with aluminium counters 
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Fig. 3 


For comparison in Fig. 3 are shown the transition curves obtained in the measure- 
ments from vertical directions by Trumpy (1943) and by Jánossy and Rossi (1940). 
Curve c has been got (by us) with a telescope of brass-counters (wall thickness 1 mm). 
It is very similar to the one found by Trumpy for vertical directions in measurements 
with brass-counters of 0,4 mm wall thickness. We have also carried out measurements 
with a telescope of aluminium-counters (of 0,6 mm wall thickness). The transition 
curve arrived at in these measurements is consistent with the curve obtained by 
Jánossy and Rossi for vertical directions (with glass counters). The observed effect 
of the shifting of the maxima of both curves is probably connected with the difference 
in the threshold energies for registering electrons in both telescopes. 

То evaluate the energies of the horizontal photons we have measured their absorp- 
tion in absorber P with a constant (1 cm thick) lead absorber S and with operating 
anticoincidence counters 4 and 5. From the examination of the absorption curve one - 
might conclude that the composition of the horizontal radiation is similar to that of 


the vertical one. 


212 M. Kogut and R. О1е5 
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To estimate the approximate intensity of the photons registered by our apparatus, 
we have used the tables of Jánossy and Messel (1950). From the position of the max- 
imum of the transition curve at X,, — 1.5 cascade units, we have estimated the average 
number Q of electrons produced by a primary photon at the maximum of the curve. 
Assuming a Poisson distribution of the actual number of electrons, we could roughly 
calculate the probability 2 of recording a photon by the detector at a depth X,, of 
the absorber, getting Q = 0.7 and A = 1 — e? = 0.5. This value is rather overesti- 
mated because of the geometrical inefficiency of the telescope. The rate of photons 
going through our telescope with aluminium-counters was about 6 per hour as comp- 
ared with the rate of ionizing particles of about 20 per hour. Our measurements of the 
absorption of the horizontal cosmic radiation show that the ratio of the number of 
mesons to the number of electrons is ~ 1, If we take into account the overestimation 
of 2, we get 1 as a roughly approximate value for the ratio of photons to electrons 
with energies exceeding some 108 eV. 

Thus in our measurements of the transition curve registering by means of pair pro- 
duction Q photons with energies of the order of 108 eV, we have not observed the 
excess of photons over electrons found by Clay for smaller energies. Besides, using 
a telescope thoroughly protected by anticoincidence-counters, we could measure, 
contrary to the statement of Rogozinski, the normal absorption of the horizontal 
cosmic radiation. Finally, the transition curve of the non-ionizing radiation registered 
in horizontal directions in our anticoincidence set was identical within the experimental 
errors with the transition curve of the photons observed in the vertical direction. 

The authors wish to thank Professor M. Miesowicz for suggesting this problem 
and for his continued interest throughout the course of this work. 
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The problem of separating the plasma motion into the so called collective part 
and the individual one, was recently the subject of several papers (Bohm and Pines 
1951, 1952, 1953, Hubbard 1954). An important step in this direction is, as it seems, 
the work of Zubariew (1953). Its fundamental idea consists in introducing into the 
wave function of the system additional variables ek, being Fourier components of the 
density operator of the system. The wave function depends on the components бк 
with indices for which k = ІЗ < k,- This cut-off procedure permits to separate out 
of the plasma motion their collective part. Parameter ky is bound in the classical theory 
with the Debye length. After such an introduction of additional variables, the Hamil- 
tonian of the system can be separated into three parts, which correspond successively to: 
free fermions, collective vibrations having a boson character and mutual interactions. 
In the event of Coulomb interaction, the second part of the Hamiltonian терелә 
oscillators vibrating with the Langmuir frequency. 

А generalization of the method of additional variables for the case of a system of 
interacting particles having different masses and charges (electrons and ions) is summa- 
ried below, and at the same time some deficiencies of Zubariew's method are removed 
thanks to an unpublished work of Migdat and Galicki, on which the present work is 
based. Our generalization permits, among others, to compute the binding energy of 
metals. The results of calculations made in the first approximation for Li, Na, K are 


given in Table I. 
The Hamiltonian of the system of interacting electrons and ions has the form: 


N, 
pi Evi ih v (|n.- ^|) - 22:3, và 


ізгі» 
ӘРІГЕ >, 7) (0) 


hPl 


(313) 
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In this formula m is the mass of the electron, r— electron coordinates, № — num- 
ber of electrons in the volume v; М, В), N,-corresponding quantities for ions. The 
wave function depends on the variables ғ, с, В, ХУ, as well as on the additional variables 
Окк, Which have to describe the collective motion of the system. Besides, the wave 
function must fulfil the additional condition 

Qix Ф (ть ©... FN, ON, R1» 2015. Во Жы t Орк) = 0 (2) 

We transform the Schrödinger equation 159/21 = Не as well as condition (2) 
by means of the unitary transformation. 
i 


y къ) 
5 = е K(K<hy) 


Р кк Ok,K 


where 


1 М N: 


ex = Үл ды ды 


N, | 9 
Око = ux e-ikrj; Р, к = — ih 
y s м. 


Instead of (2) we obtain the condition 


(бк — 0509 = 0 (4) 
We separate one part of the Hamiltonian, corresponding to the Coulomb interac- 
tion into a *far^ and а “near“ interaction (k X k,, k > kj), for instance 


e (kej ERp ЕР (3) 


in particular 


1 
z У и (= о Oe 2 22 (9) + 


2 44 № K(k < ke) k (k <ko) 
2 „їй, - гу) ©) 
Аа к>) 29 


After performing this transformation we introduce new variables by 


А 7 4 
Оһо = V22n0 qro; Pko = L——— 70; Ako = GNI 
V2 Axo 4т М, v(k) 
5 F . 
Qox = У gox; Pox = —-— Лок; An = : NG MV (6) 
V2 Ave AMN; (E) 


Terms corresponding to the collective vibrations (k <k,), after making use of 
condition (4), can be expressed in the form 


hw hw 
Н,, == ES (%x,0%—x,;0 + 9к,о q—x,0) + TOS (лок 70,-к-- qok qo,-X) + 


k(R< ko) 
— (ko 9,-к + 9-0 ды | i 
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Hence it сап be seen, that the collective vibrations correspond with vibrations of 
: DES. 
coupled oscillators, with the coupling constant A = > (Фи, where соу = 


Anne? Алпе? : ; 
= ----, 09 = — —- are Langmuir frequencies of electrons and ions, 
m M 
and лу, ng — densities of electrons and ions. 
From (7) we may find the | binding energy of the ground state. For this purpose we 


pass to normal coordinates л, q, in which the system (7) can be represented as a E 


of not interacting. oscillators, vibrating at frequencies: о’ = Vo? + о: + 02, 0" = 0, 
£ ID M CRT ze 
Нос = X > (x лк + qr q-k) (8) 
k (k< ko) 


After passing to the boson occupation-number representation, we have: 


Н ms ho’ P3 (bg by, sp У;) (9) 
k(k< ko) 


Hence we get the energy of the ground state 
En ocn (10) 


where п’ is the number of independent additional coordinates. We assume that n’~ 


N, = N, (Bohm and Pines 1953). If in the Hamiltonian (7) we omit the coupling, then 
the energy of the basic state will be E n'h (оу + а»). We obtain the binding energy 
W from the relation: 

n'h 


т (о, + о) — W = eve (11) 


Detailed consideration and examination of the remaining terms of the Hamilto- 
nian wil be given in the main paper. 


Table I. Energy W in kcal/mol 


Experimental W calculated from W ш Focks W in Fock's 
W (11) approximation approximation, 
taking into the 
correlation 
Li 39 28 14,6 | -86,2 
Ха 26 10 ЕЛ! 24,5 
EK 22,6 6,5 —0,7 16,5 
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